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Abstract 

In  collaboration  with  AFRL  colleagues  we  have  developed  a  methodology  for  measuring  localized 
microstructure  metrics  in  dendritic  microstructures  such  as  single  crystal  nickel-based  superalloys, 
with  potential  application  to  a  wide  range  of  material  systems.  The  primary  dendrite  arm  spacing 
is  correlated  to  processing  (solidification  rate),  the  micro  structure  of  the  material  (interdentritic 
eutectic  particles  and  voids),  and  properties  (fatigue  behavior  and  creep  strength).  In  this  project 
a  Voronoi-based  approach  for  spatial  point  pattern  analysis  has  been  applied  to  experimental  and 
synthetic  dendritic  microstructures.  This  technique  was  used  to  quantify  the  distribution  of  local 
primary  dendrite  arm  spacings,  their  spatial  distribution,  and  their  correlation  with  interdentritic 
eutectic  particles.  Several  peer-reviewed  journal  articles,  including  an  Editor’s  Choice  article 
published  in  Metallurgical  and  Materials  Transactions,  and  conference  presentations  describe  this 
in  more  detail.  Furthermore,  to  transfer  this  technology,  a  GUI  based  computational  too  has  been 
developed  to  facilitate  the  use  of  these  metrics  in  synthetic  and  experimental  microstructures,  with 
eventual  application  to  link  processing  with  structural  properties  based  on  the  complex  turbine  blade 
microstructure. 
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1  Final  Project  Summary 

The  following  section  includes  details  of  the  proposal  and  program  plans  for  the  AFOSR  project. 

1.1  Research  Effort 

This  3-year  research  project  qualitatively  and  quantitatively  mapped  the  relationship  between  the 
hierarchical  microstructure  and  the  elevated  temperature  creep  deformation  and  failure  behavior  of 
single  crystal  nickel-based  superalloys.  This  research  established  a  methodology  for  incorporating 
microstructure-based  structure-property  relationships  into  mechanics-based  material  models  and 
answering  questions  related  to  “how  microstructural  features  alter  localized  deformation  at  elevated 
temperatures”  and  “what  is  the  materials  uncertainty  associated  with  elevated  temperature  creep 
deformation”  for  developing  the  next  generation  of  elevated  temperature  structural  materials  for 
propulsion  systems  applications  through  materials-based  design  approaches. 

On  a  large  scale,  a  problem  central  to  designing  the  next  generation  propulsion  systems  for 
the  USAF  is  that  the  present  length  of  time  required  for  new  material  design  is  much  larger  than 
the  timelines  required  for  development  of  commercial  engines  and  engine  components.  At  the 
Superalloys  2008  conference,  the  keynote  address  postulated  that  the  grand  challenge  for  materials 
development  is  to  drastically  reduce  this  development  time  for  new  materials  and  processes  without 
adding  development  risk.  Indeed,  accelerating  the  pace  of  discovery,  development,  manufacturing, 
and  deployment  of  advanced  materials  systems  has  been  described  as  crucial  to  achieving  global 
competitiveness  in  the  21st  century  in  a  “Materials  Genome  Initiative”  white  paper  recently  released 
by  the  White  House’s  National  Science  and  Technology  Council.  Moreover,  recent  discussions  and 
reports  on  integrated  computational  materials  engineering  (ICME)  provide  a  conceptual  roadmap 
for  how  to  design  materials,  but  the  community  also  agrees  that  there  are  still  a  number  of  challenges 
that  must  be  addressed  before  ICME  reaches  its  full  potential.  Nevertheless,  success  stories  in  the 
automotive  and  aerospace  sectors  show  that  the  aforementioned  grand  challenge  is  solvable. 

Accelerating  the  pace  of  discovery  requires  physically -based  models  that  describe  the  “processing- 
structure-property-performance”  relationships.  As  Figure  1  shows,  these  models  require  information 
about  how  material  processing/chemistry  relates  to  structure,  how  structure  relates  to  properties, 
and  how  properties  ultimately  relate  to  the  performance  of  a  specific  component  for  a  particular 
application.  As  an  example,  the  /  precipitate  micro  structure  (structure,  left)  is  shown  in  Figure 
1  with  results  from  elevated  temperature  creep  testing  (property,  right).  The  structure-property 
relationship  (middle)  captures  the  link  between  the  volume  fraction  of  /  precipitates  and  creep 
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rupture  life.  Quantifying  these  relationships  is  invaluable  both  for  designing  new  materials  and  for 
optimizing  the  properties/performance  of  existing  materials.  Developing  the  models,  databases,  and 
tools  for  developing  these  relationships  will  help  establish  the  infrastructure  required  to  solve  this 
complex  engineering  problem. 

On  a  smaller  scale,  a  relevant  problem  for  the  USAF  is  that  increasingly  higher  temperatures 
and  pressures  in  the  high  pressure  turbine  section  of  the  modern  turbine  engine  leads  to  a  steep 
decrease  in  properties  at  temperatures  above  850  °C  for  single  crystal  nickel-based  superalloys. 
However,  to  solve  this  specific  problem  requires  developing  mechanics-based  material  models  with 
uncertainty  that  can  capture  the  appropriate  structure-property  relationships  by:  (1)  measuring 
details  of  the  hierarchical  microstructure  and  its  evolution  and  (2)  relating  this  to  properties  at 
elevated  temperatures  under  extreme  conditions  (i.e.,  creep,  low-  and  high-cycle  fatigue,  fretting, 
embrittlement,  and  environment  effects).  In  this  sense,  these  models  will  relate  the  initial  material 
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Figure  1.  Schematic  of  the  “processing-structure -property -performance”  materials  design  concept.  The  pro¬ 
posed  research  will  focus  on  structure-property  relationships  in  single  crystal  Ni-based  superalloys 
for  materials-based  design.  The  example  shows  how  precipitate  volume  fraction  affects  creep 
rupture  life. 
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microstructure  to  the  properties  that  ultimately  influence  performance  of  components  and  systems. 
To  do  so  requires  that  the  relevant  features  and  their  variability  are  captured  via  statistical  descriptors. 
Moreover,  the  relationship  between  structure  and  properties  must  be  firmly  established  and  carefully 
validated  with  experiments  in  order  to  not  increase  development  risk.  There  are  multiple  challenges 
to  developing  microstructure -based  predictive  models  and  intrinsically  coupling  with  variability 
associated  with  microstructure  heterogeneities: 

1.  What  metrics  should  be  used  to  describe  a  digital  representation  of  material  micro  structure? 
How  does  airfoil  geometry  (processing)  affect  microstructure?  What  is  the  uncertainty 
associated  with  microstructure  statistics  (i.e.,  representative  volume  size)?  How  does  material 
micro  structure  evolve  during  service? 

2.  How  does  micro  structure  feature  heterogeneity  relate  to  local  deformation  character?  Can 
in-situ  experiments  be  used  to  provide  local  structure-property  relationships?  How  does  the 
microstructure  affect  cumulative  damage  during  elevated  temperature  deformation? 

3.  What  is  the  link  between  extreme  values  of  creep  rupture  life  and  micro  structure?  What  is 
the  materials  uncertainty  associated  with  creep  properties?  Can  this  be  related  to  the  local 
microstructure? 

To  address  these  challenges,  this  research  uses  a  combination  of  microscopy  and  elevated 
temperature  mechanical  testing,  combined  with  computational  approaches.  Here,  a  unique,  multi¬ 
pronged  approach  is  pursued  to  understand  the  role  of  material  micro  structure  and  their  associated 
uncertainties  on  the  deformation  behavior  of  aerospace  materials  at  a  number  of  length  scales 
that  span  from  the  macroscale  (homogeneous  microstructure  and  properties)  to  the  microscale 
(containing  resolution  of  individual  microstructure  features).  In  particular,  this  research  uses  Air 
Force-relevant  single  crystal  nickel-based  superalloys  to  explore  the  influence  of  micro  structure  on 
creep  deformation. 

1.2  Significance  and  Relevance  to  USAF 

Single  crystal  nickel  turbine  blades  are  being  utilized  in  turbine  jet  engines  throughout  industry 
because  of  their  superior  creep,  stress  rupture,  melt  resistance,  and  thermomechanical  fatigue 
capabilities  over  poly  crystalline  alloys.  With  respect  to  its  significance  to  the  USAF,  single  crystal 
alloys  are  used  for  turbine  blades  within  the  high  pressure  turbine  (HPT)  section  of  turbine  engines, 
e.g.,  in  the  GE  FI  10  engine.  The  FI  10  turbine  engine  has  been  the  engine  of  choice  for  the  F-16 
since  there  was  an  engine  choice  for  that  aircraft;  fully  86%  of  the  USAF  F-16C/Ds  and  75%  of  all 
front  line,  combat-coded  F-16s  are  powered  by  the  GE  FI  10.  Part  of  the  motivation  for  choosing 
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this  particular  alloy  class  is  that  the  limiting  factor  for  specific  thrust  within  the  turbine  engine 
is  the  temperature  of  the  turbine  section  (where  the  single  crystal  turbine  blades  reside).  Thus,  a 
fundamental  understanding  of  the  hierarchical  microstructure  influences  that  contribute  to  failure 
can  contribute  to  both  an  improvement  in  the  ability  to  predict  life  in  these  systems  as  well  as  the 
potential  to  re-engineer  these  alloys  or  redesign  the  components  for  improved  elevated  temperature 
properties  and  propulsion  system  efficiency. 

As  the  operating  temperatures  in  the  hot  section  of  the  turbine  engine  continue  to  increase,  the 
airfoil  design  and  geometry  is  becoming  increasingly  complex  -  wall  thicknesses  are  decreasing  and 
cooling  passages  are  more  intricate  to  keep  the  turbine  blade  cool.  Therefore,  it  has  become  even 
more  vital  to  understand  how  microstructure  affects  high  temperature  properties  as  wall  thicknesses 
and  geometry  within  the  airfoil  approach  the  length  scales  of  some  of  the  microstructure  features 
themselves.  This  proposal  will  utilize  experiments  for  testing  specimens  that  are  on  the  order  of  the 
wall  thickness  of  turbine  blades,  as  it  has  previously  been  found  that  specimen  thickness  plays  a 
commanding  role  in  elevated  temperature  creep  properties.  Creep  deformation  testing  is  especially 
relevant  because  the  severe  loads  experienced  in  long  cruise  missions  have  a  large  fraction  of  hot 
dwell  time.  Again,  this  testing  is  motivated  by  mission  conditions  experienced  by  thin  sections 
within  HPT  blades. 

1.3  Final  Report  Summary  of  Accomplishments 

The  summary  and  accomplishments  of  the  current  AFOSR  project  is  as  follows: 

1.  In  collaboration  with  AFRL  colleagues,  we  have  developed  a  methodology  for  measuring 
localized  microstructure  metrics  in  dendritic  microstructures  (e.g.,  single  crystal  nickel -based 
superalloys),  with  potential  application  to  order/disorder  in  a  wide  range  of  material  systems. 
Characterizing  the  primary  dendrite  arm  spacing  in  directionally- solidified  microstructures 
is  an  important  step  for  developing  process-structure-property  relationships  by  enabling  the 
quantification  of  (i)  the  influence  of  processing  on  microstructure  and  (ii)  the  influence  of 
microstructure  on  properties.  Thin-walled  directionally- solidified  structures  (e.g.,  a  turbine 
blade)  require  new  approaches  for  characterizing  the  dendrite  arm  spacing  and  the  microstruc¬ 
ture.  In  this  work,  we  utilized  a  new  Voronoi-based  approach  for  spatial  point  pattern  analysis 
that  was  applied  to  an  experimental  dendritic  micro  structure.  This  technique  utilizes  a  Voronoi 
tessellation  of  space  surrounding  the  dendrite  cores  to  determine  nearest  neighbors  and  the 
local  primary  dendrite  arm  spacing.  In  addition,  we  compared  this  technique  to  a  recent 
distance-based  technique,  the  Warnken-Reed  method,  and  a  modification  to  this  using  Voronoi 
tessellations,  along  with  the  minimal  spanning  tree  method.  Moreover,  a  convex  hull-based 
technique  was  used  to  include  edge  effects  for  such  techniques,  which  can  be  important  for 
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thin  specimens.  These  methods  were  used  to  quantify  the  distribution  of  local  primary  dendrite 
arm  spacings  as  well  as  their  spatial  distribution  for  an  experimental  directionally-solidified 
superalloy  micrograph.  Last,  eutectic  particles  were  segmented  to  correlate  distances  from 
dendrite  cores  and  Voronoi  vertices  to  the  occurrence  and  size  of  these  interdendritic  features. 
Interestingly,  with  respect  to  the  distance  from  the  dendrite  core,  it  was  found  that  there  is  a 
greater  probability  of  occurrence  of  large  eutectic  particles  (>  410  jUm)  over  small  particles  at 
distances  greater  than  approximately  1/3  of  the  bulk-measured  primary  dendrite  arm  spacing. 
This  systematic  study  of  the  different  techniques  for  quantifying  local  primary  dendrite  arm 
spacings,  and  their  effect  on  microstructure,  can  be  an  important  step  for  correlating  with 
both  processing  and  properties  in  single  crystal  nickel-based  superalloys. 

2.  Furthermore,  to  transfer  this  technology,  a  GUI-based  computational  tool  has  been  developed 
to  facilitate  the  use  of  these  metrics  in  synthetic  and  experimental  microstructures,  with 
eventual  application  to  link  processing  with  structural  properties  based  on  the  complex 
turbine  blade  micro  structure.  This  tool  is  being  used  to  examine  the  uncertainty  with  various 
methods  for  measuring  local  micro  structure  characteristics  and  to  examine  the  sensitivity  of 
using  such  techniques  for  evaluating  microstructure.  In  addition  to  the  GUI  tool,  a  manual 
describing  its  use  has  been  produced  and  will  be  released  along  with  the  tool  after  publication 
of  the  first  combined  study  with  synthetic  and  experimental  microstructures  (anticipated  Fall 
2014).  Additionally,  further  techniques  for  characterizing  micro  structure,  such  as  the  minimal 
spanning  tree  methodology,  were  explored  in  collaboration  with  Kiran  Solanki  at  Arizona 
State  University. 

3.  At  the  Army  Research  Laboratory  (ARL)  where  the  original  PI  for  this  project,  Mark  A. 
Tschopp,  is  a  researcher  now,  a  small  scale  (down  to  150-/im  diameter)  high  temperature 
(<1200  °C)  shear  punch  tester  was  built  to  examine  the  utility  of  rapid,  small- volume,  high 
temperature  creep  uncertainty  studies  on  aerospace  alloys.  As  most  creep  tests  can  take  several 
days  on  bulk  specimens,  this  may  allow  for  more  rapid  capability  for  assessing  Larson-Miller 
plots  as  well  as  providing  the  potential  capability  to  examine  the  creep  response  in  very  thin 
sections  such  as  in  the  actual  turbine  blade. 

4.  Mark  A.  Tschopp  was  the  lead  organizer  for  a  symposium  on  “Multi-scale  Microstructure, 
Mechanics  &  Prognosis  of  High  Temperature  Alloys”  for  TMS  2015  (March  15-19,  2015  in 
Orlando,  Florida)  with  co-organizers  Jeffrey  Evans,  Jonathan  Cormier,  and  Qiang  (Charles) 
Feng.  This  symposium  was  sponsored  by  the  High  Temperature  Alloys  Committee  and 
was  co-sponsored  by  the  Computational  Materials  Science  committee  and  the  Mechanical 
Behavior  of  Materials  committee.  There  were  62  presentations  in  this  symposium,  including 
one  from  authors  on  this  project. 
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5.  Some  of  the  results  of  this  work  have  been  featured  in  an  Editor’s  Choice  article  published  in 
Metallurgical  and  Materials  Transactions  as  well  as  multiple  talks,  most  notably  invited  talks 
at  the  annual  TMS  2013  and  ASME  IMECE  2013  conferences. 

6.  Additionally,  recent  discussions  with  the  Propulsion  branch  at  the  Army  Research  Laboratory 
may  help  in  establishing  future  opportunities  for  collaboration  between  the  Army  and  the  Air 
Force  in  the  area  of  high  temperature  aerospace  materials  research. 

1.4  Papers  and  Presentations  delivered 

Peer  Reviewed,  Archived  Journal  Articles 

1.  M.A.  Tschopp,  A.L.  Oppedal,  J.D.  Miller,  M.A.  Groeber,  A.H.  Rosenberger,  K.N.  Solanki, 
Characterizing  primary  dendritic  microstructures  to  quantify  the  processing-structure -property 
relationship  in  single  crystal  nickel-based  superalloys,  ed.  J.  Y.  Hwang,  C.  Bai,  J.  Carpenter, 
S.  Ikhmayies,  B.  Li,  N.  Monteiro,  Z.  Peng,  M.  Zhang,  Characterization  of  Minerals,  Metals, 
and  Materials  2013,  TMS  (The  Minerals,  Metals,  &  Materials  Society),  2013:  pp.  301-310. 

2.  This  paper  was  selected  as  an  Editor’s  Choice  as  having  special  significance. 

M.A.  Tschopp,  J.D.  Miller,  A.L.  Oppedal,  K.N.  Solanki,  Characterizing  the  Local  Primary 
Dendrite  Arm  Spacing  in  Directionally  Solidified  Dendritic  Microstructures,  Metallurgical 
and  Materials  Transactions  A.  45  (2014)  426-437. 

3.  M.A.  Tschopp,  J.D.  Miller,  A.L.  Oppedal,  K.N.  Solanki,  Evaluating  Local  Primary  Dendrite 
Arm  Spacing  Characterization  Techniques  Using  Synthetic  Directionally  Solidified  Dendritic 
Microstructures,  Metallurgical  and  Materials  Transactions  A.  46  (2015)  4610-4628. 

Peer  Reviewed  Abstract,  Invited  Presentation 

1.  M.A.  Tschopp,  A.L.  Oppedal,  J.D.  Miller,  M.A.  Groeber,  A.H.  Rosenberger,  K.N.  Solanki, 
Characterizing  primary  dendritic  microstructures  to  quantify  the  processing-structure -property 
relationship  in  single  crystal  nickel-based  superalloys,  in:  TMS  2013  Annual  Meeting,  March 
3-7,  2013,  San  Antonio,  TX 

2.  M.A.  Tschopp,  A.L.  Oppedal,  J.D.  Miller,  M.A.  Groeber,  A.H.  Rosenberger,  K.N.  Solanki,  Un¬ 
derstanding  the  Processing-Structure-Property  Relationships  in  Single  Crystal  Nickel-Based 
Superalloys,  in  ASME  2013  International  Mechanical  Engineering  Congress,  November 
15-21,  2013,  San  Diego,  CA 
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3.  M.A.  Tschopp,  A.L.  Oppedal,  K.N.  Solanki,  Local/Global  Measurement  of  Primary  Dendrite 
Arm  Spacing  in  Single  Crystal  Nickel-Based  Superalloys, in:  TMS  2015  Annual  Meeting, 
March  15-19,  2013,  Orlando,  FL 

Program  Review  Meetings 

1.  M.A.  Tschopp,  A.L.  Oppedal,  S.  Turnage,  (Poster)  Hierarchically-driven  approach  for  quanti¬ 
fying  materials  uncertainty  in  creep  deformation  and  failure  of  aerspace  materials,  Multi-scale 
Structural  Mechanics  and  Prognosis  Awardee  Review  Meeting,  July  25,  2012,  Washington, 
DC 

2.  M.A.  Tschopp,  A.L.  Oppedal,  R.  Carino,  S.  Turnage,  K.  N.  Solanki,  (Presentation)  Hierarchically 
driven  approach  for  quantifying  materials  uncertainty  in  creep  deformation  and  failure  of 
aerspace  materials,  Multi-scale  Structural  Mechanics  and  Prognosis  Awardee  Review  Meeting, 
July  24,  2013,  Washington,  DC 

3.  M.A.  Tschopp,  A.L.  Oppedal,  R.  Carino,  S.  Turnage,  K.  N.  Solanki,  (Presentation)  Hierarchically 
driven  approach  for  quantifying  materials  uncertainty  in  creep  deformation  and  failure  of 
aerspace  materials,  Multi-scale  Structural  Mechanics  and  Prognosis  Awardee  Review  Meeting, 
September  5,  2014,  Albuquerque,  NM 
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Characterizing  the  Local  Primary  Dendrite  Arm  Spacing 
in  Directionally  Solidified  Dendritic  Microstructures 


MARK  A.  TSCHOPP,  JON  D.  MILLER,  ANDREW  L.  OPPEDAL, 
and  K1RAN  N.  SOLANKI 

Characterizing  the  spacing  of  primary  dendrite  arms  in  directionally  solidified  microstructures  is 
an  important  step  for  developing  process-structure-property  relationships  by  enabling  the 
quantification  of  (i)  the  influence  of  processing  on  microstructure  and  (ii)  the  influence  of 
microstructure  on  properties.  In  this  work,  we  utilized  a  new  Voronoi-based  approach  for 
spatial  point  pattern  analysis  that  was  applied  to  an  experimental  dendritic  microstructure.  This 
technique  utilizes  a  Voronoi  tessellation  of  space  surrounding  the  dendrite  cores  to  determine 
nearest  neighbors  and  the  local  primary  dendrite  arm  spacing.  In  addition,  we  compared  this 
technique  to  a  recent  distance-based  technique  and  a  modification  to  this  using  Voronoi  tes¬ 
sellations.  Moreover,  a  convex  hull-based  technique  was  used  to  include  edge  effects  for  such 
techniques,  which  can  be  important  for  thin  specimens.  These  methods  were  used  to  quantify  the 
distribution  of  local  primary  dendrite  arm  spacings,  their  spatial  distribution,  and  their  corre¬ 
lation  with  interdendritic  eutectic  particles  for  an  experimental  directionally  solidified  Ni-based 
superalloy  micrograph.  This  can  be  an  important  step  for  correlating  processing  and  properties 
in  directionally  solidified  dendritic  microstructures. 

DOL  10.1007/sl  1661-01 3-1985-3 
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I.  INTRODUCTION 

Developing  an  enhanced  understanding  of 
mechanical  behavior  in  materials  relies  upon  sufficiently 
characterizing  microstructure  details  at  the  relevant 
length  scales  that  contribute  to  this  behavior.  Moreover, 
to  truly  enhance  the  predictive  capability  of  processing- 
structure-property  models  that  aim  to  improve  material 
performance  requires  a  quantitative  stereological 
description  of  the  relevant  microstructure  features  and, 
thereby,  the  material  itself.  Predictive  models  that 
effectively  capture  the  linkage  between  processing  and 
properties  (through  microstructure)  can  be  utilized 
within  an  integrated  computational  materials  engineer¬ 
ing  (ICME)  approach  to  design  materials  and  accelerate 
their  insertion  into  application. 

The  focus  of  the  present  work  is  on  single-crystal 
nickel-based  superalloys,  which  are  used  in  turbine 
blades  within  the  high  temperature  section  of  the 
modern  turbine  engine!11,21  In  single-crystal  nickel-based 
superalloys,  there  are  a  number  of  length  scales  of 
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microstructure  that  contribute  to  mechanical  behavior, 
ranging  from  the  /  precipitates  to  pores  and  eutectic 
particles  to  the  dendrites  themselves.  At  the  largest 
microstructure  length  scale  in  directionally  solidified 
single-crystal  microstructures,  the  features  of  interest  are 
the  dendrites;  many  features  at  lower  length  scales  (e.g., 
eutectic  particles,  precipitates,  etc.)  or  at  similar  scales 
(e.g.,  porosity,  freckle  defects,  etc.)  are  strongly  associated 
with  the  dendrite  arm  spacing  and  morphology.13  71  The 
solidification  morphology  associated  with  dendrite  arm 
spacing  has  been  described  in  the  literature.18,91  Histori¬ 
cally,  the  primary  dendrite  arm  spacing  (PDAS)  has  been 
found  to  correlate  with  processing  (e.g.,  solidification 
rate)17,10  141  as  well  as  with  properties  (e.g.,  creep  strength, 
fatigue  properties). 115,61  For  instance,  Lamm  and  Singer161 
produced  single-crystal  nickel-based  microstructures 
(PWA  1483)  with  a  varied  range  of  different  dendrite 
arm  spacings  (250  to  600  fim)  and  found  that  decreasing 
the  mean  dendrite  arm  spacing  was  associated  with  an 
increased  high-cycle  fatigue  life.  The  fatigue  cracks  were 
found  to  originate  at  shrinkage  porosity  and  the  largest 
pores  correlated  with  a  large  PDAS. 

The  traditional  approach  for  measuring  PDAS  in 
single-crystal  metals,  whereby  the  number  of  dendrite 
cores  in  a  specified  area  is  related  to  the  dendrite  arm 
spacing1 10,1 6,1 71  is  given  by: 


where  X  is  PDAS,  A  is  the  area  analyzed,  n  is  the  number 
of  dendrites,  and  c  is  a  coefficient  that  depends  on  the 
microstructure.  McCartney  and  Hunt1101  showed  that 
c  =  0.5  for  a  random  array  of  points,  c  =  1  for  a 

METALLURGICAL  AND  MATERIALS  TRANSACTIONS  A 


Figure  2.  Example  peer-reviewed  journal  publication  published  in  the  current  AFOSR  project. 
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Figure  3.  Example  invited  presentation  on  materials  characterization  for  superalloys  delivered  at  national  and 
international  professional  society  conferences. 
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Figure  4.  Example  Multi-scale  Structural  Mechanics  and  Prognosis  Awardee  Review  presentation  . 
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Meeting:  2015  TMS  Annual  Meeting  &  Exhibition 
Symposium: 

Multi-scale  Microstructure,  Mechanics  &  Prognosis  of  High  Temperature  Alloys 
Sponsorship: 

High  Temperature  Alloys  Committee  (Sponsor) 

Mechanical  Behavior  (co-Sponsor),  Computational  Materials  Science  (co-Sponsor) 

Organizers: 

Mark  A.  Tschopp,  U.S.  Army  Research  Laboratory 
Jeffrey  L.  Evans,  University  of  Alabama  in  Huntsville 

Jonathan  Cormier,  Ecole  Nationale  Superieure  de  Mecanique  et  d'Aerotechnique  (ENSMA),  France 
Qiang  (Charles)  Feng,  University  of  Science  and  Technology  Beijing  (USTB),  China 

Scope: 

A  number  of  critical  applications  in  industry  demand  the  use  of  high  temperature  alloys  that  can 
withstand  various  extreme  environments  under  elevated  temperature  conditions.  Reliably  predicting  the 
life  of  these  components,  which  may  be  subjected  to  damage  associated  with  creep  deformation,  cyclic 
loading,  environmental  degradation,  and  combinations  thereof,  is  a  challenging  endeavor.  Under  these 
conditions,  the  state  of  the  microstructure  after  processing  and  its  evolution  during  service  lead  to 
damage  mechanisms  that  manifest  over  multiple  length  scales,  ranging  from  quantum  and  atomistic 
scales  to  the  mesoscale  to  the  macroscale. 

This  symposium  will  provide  a  venue  for  presenting  recent  achievements  in  understanding  the 
microstructure  evolution  and  mechanical  behavior  in  high  temperature  alloys  over  multiple  scales  in 
order  to  ultimately  predict  the  prognosis  and  life  of  components.  Advances  in  experimental  and 
computational  capability  have  greatly  improved  our  ability  to  understand  and  quantify  deformation 
mechanisms  in  high  temperature  alloys.  In  particular,  a  goal  of  this  symposium  is  to  accelerate  the 
development  and  acceptance  of  new  methodologies  for  improving  prognosis  through  understanding  the 
fundamental  relationships  between  material  microstructure  and  mechanical  behavior  in  these  alloys.  It 
is  expected  that  this  symposium  will  include  talks  ranging  from  atomistic,  discrete  dislocation,  and 
continuum  mechanics  approaches  for  various  length  scales  as  well  as  experimental  mechanics  results 
that  elucidate  the  behavior  of  these  alloys. 

The  subject  areas  of  this  symposium  include,  but  are  not  limited  to: 

•  Multi-scale  modeling  of  high  temperature  deformation  and  damage 

•  Single  crystal  and  polycrystal  plasticity  models 

•  Deformation  and  damage  based  life  prediction  techniques 

•  Modeling  and  experimental  approaches  to  creep-fatigue-environment  interactions 

•  Discrete  dislocation  dynamics  and  mesoscale  (phase  field)  modeling  of  creep  defonnation 

•  Experimental  methods  for  materials  prognosis  and  structural  health  monitoring 

•  Understanding  and  quantifying  microstructure-property  relationships  at  high  temperature 

•  Computational  and  experimental  approaches  for  accelerated  materials  design  of  high  temperature 
alloys 

•  Effect  of  chemistry  and  processing  on  high  temperature  structure  and  properties 


Figure  5.  Symposium  on  high  temperature  microstructure  and  mechanics  of  aerospace  materials  organized 
for  TMS  2015. 
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2  TMS  Conference  Paper,  March  2013 


M. A.  Tschopp,  A.L.  Oppedal,  J.D.  Miller,  M.A.  Groeber,  A.H.  Rosenberger,  K.N.  Solanki,  Charac¬ 
terizing  primary  dendritic  microstructures  to  quantify  the  processing-structure-property  relationship 
in  single  crystal  nickel-based  superalloys,  ed.  J.  Y.  Hwang,  C.  Bai,  J.  Carpenter,  S.  Ikhmayies,  B.  Li, 

N.  Monteiro,  Z.  Peng,  M.  Zhang,  Characterization  of  Minerals,  Metals,  and  Materials  2013,  TMS 
(The  Minerals,  Metals,  &  Materials  Society),  2013:  pp.  301-310. 
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Abstract 

Characterizing  the  spacing  of  primary  dendrite  arms  in  directionally-solidified  microstructures  is 
an  important  step  for  developing  process-structure-property  relationships  by  enabling  the 
quantification  of  (i)  the  influence  of  processing  on  microstructure  and  (ii)  the  influence  of 
microstructure  on  properties.  The  research  objective  herein  is  to  evaluate  the  capability  of 
various  conventional  approaches,  as  well  as  new  or  modified  approaches,  for  spatial  point  pattern 
analysis  with  application  to  characterizing  experimental  dendritic  microstructures.  Both 
computer-generated  and  experimental  dendritic  microstructures  are  used  for  this  analysis  along 
with  numerous  techniques  based  on  the  nearest  neighbor  spacing,  Voronoi  tessellation,  Delaunay 
triangulation,  or  graph  theory.  Comparison  of  new  metrics  with  traditional  primary  dendrite  arm 
spacing  metrics  will  also  be  discussed  for  both  local  and  global  measures.  The  current  methods 
investigated  will  supply  information  of  local  spacing  and  coordination  number  while  addressing 
edge  effects,  parameter  sensitivity,  and  correlation  with  interdendritic  features,  thus  providing 
insight  into  how  processing  affects  properties. 

Introduction 

Developing  an  enhanced  understanding  of  mechanical  behavior  in  materials  relies  upon 
sufficiently  characterizing  microstructure  details  at  the  relevant  length  scales  that  contribute  to 
this  behavior.  Moreover,  to  truly  enhance  the  predictive  capability  of  processing-structure- 
property  models  that  aim  to  improve  material  performance  requires  a  quantitative  description  of 
the  relevant  microstructure  features  and,  thereby,  the  material  itself. 
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In  this  work,  we  focus  on  single  crystal  nickel-based  superalloys,  which  are  used  in  turbine 
blades  within  the  high  temperature  section  of  the  modem  turbine  engine  (Reed,  2006  ).  In  single 
crystal  nickel-based  superalloys,  there  are  a  number  of  length  scales  of  microstructure  that 
contribute  to  mechanical  behavior,  ranging  from  the  primary,  secondary,  and  tertiary  gamma 
prime  precipitates  to  pores  and  eutectic  particles  to  the  dendrites  themselves.  As  the  dendritic 
microstructure  lies  at  the  largest  length  scales,  it  is  important  to  characterize  their  spacing,  as 
historically  the  PDAS  (primary  dendrite  arm  spacing)  has  been  found  to  correlate  with 
processing  (e.g.,  solidification  rate)  (Wang  et  al.,  2003;  Brundidge,  2011)  as  well  as  with 
properties  (e.g.,  creep  strength,  fatigue  properties)(Wilson  et  al.,  2008;  Lamm  and  Singer,  2007). 
Characterizing  the  spacing  of  primary  dendrite  arms  in  directionally-solidified  microstructures  is 
an  important  step  for  developing  process-structure-property  relationships  by  enabling  the 
quantification  of  (i)  the  influence  of  processing  on  microstructure  and  (ii)  the  influence  of 
microstructure  on  properties. 

The  traditional  approach  for  measuring  primary  dendrite  arm  spacing  in  single  crystal  metals, 
whereby  the  number  of  dendrite  cores  in  a  specified  area  is  related  to  the  dendrite  arm  spacing 
(Flemings,  1974;  Jacobi  and  Schwerdtfeger,  1976;  Wamken  and  Reed,  2010a,  Wamken  and 
Reed  2010b): 

where  X  is  primary  dendrite  arm  spacing,  A  is  the  area  analyzed,  n  is  the  number  of  dendrites,  and 
c  is  a  coefficient  that  depends  on  the  microstructure  -  is  insufficient  for  capturing  local  arm 
spacings  or  the  dendrite  arm  spacing  distribution,  and  may  provide  problems  with  complex 
geometries  such  as  turbine  blades. 

The  research  objective  herein  is  to  evaluate  the  capability  of  various  conventional  approaches,  as 
well  as  new  or  modified  approaches,  for  spatial  point  pattern  analysis  with  application  to 
characterizing  experimental  dendritic  microstructures.  Both  computer-generated  and 
experimental  dendritic  microstructures  are  used  for  this  analysis  along  with  numerous  techniques 
based  on  the  nearest  neighbor  spacing,  Voronoi  tessellation,  Delaunay  triangulation,  or  graph 
theory.  Comparison  of  new  metrics  with  traditional  primary  dendrite  arm  spacing  metrics  will 
also  be  discussed  for  both  local  and  global  measures.  The  current  methods  investigated  will 
supply  information  of  local  spacing  and  coordination  number  while  addressing  edge  effects, 
parameter  sensitivity,  and  correlation  with  interdendritic  features,  thus  providing  insight  into 
how  processing  affects  properties. 
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Methodology 

The  approach  utilized  in  this  work  to  measure  the  local  dendrite  arm  spacing  is  based  on  a 
Voronoi  tessellation  of  the  spatial  array  of  dendrite  cores.  To  illustrate  how  the  present  method 
works  and  differs  from  some  other  published  methods,  we  generated  a  synthetic  5x5  cubic 
pattern  of  points  with  a  small  degree  of  noise  (100%  noise  fraction,  0.20ao  noise  fraction  [2]  ),  as 
shown  in  Figure  la.  For  the  purposes  of  describing  several  different  methods  shown  in  Figure  1, 
this  synthetic  pattern  of  points  can  be  considered  as  the  cores  of  primary  dendrites. 

One  such  method  for  measuring  the  local  dendrite  arm  spacing  is  the  Wamken-Reed  method 
(Wamken  and  Reed,  2011).  The  Wamken-Reed  method  calculates  the  dendrite  arm  spacing  for 
a  single  point  (black  dot)  by  starting  with  an  initial  number  of  nearest  neighbors  (3  closest 
neighbors)  and  iteratively  adding  potential  nearest  neighbors  that  are  within  a  cutoff  distance 
defined  by  the  already-added  nearest  neighbors.  For  instance,  the  inner  circle  in  Figure  lb 
represents  the  average  spacing,  dave,  of  these  neighbors  and  the  outer  circle  represents  the  cutoff 
for  adding  the  next  neighbor,  dave  +a*dstd-  Neighbors  continue  to  be  added  until  the  cutoff  does 
not  include  any  new  neighbors.  The  local  coordination  number  and  dendrite  arm  spacing  is 
calculated  from  the  neighbors  added  (shown  as  red  dots).  Flowever,  if  the  standard  deviation  of 
the  distances  of  the  nearest  neighbors  is  large,  this  technique  can  continue  to  add  nearest 
neighbors  beyond  the  first  nearest  neighbors.  Clearly,  a  method  for  restricting  the  number  of 
nearest  neighbors  using  such  a  technique  is  necessary. 

A  simple  way  of  identifying  the  potential  first  nearest  neighbors  is  to  perform  a  Voronoi 
tessellation  of  the  space  surrounding  the  points,  as  shown  in  Figure  lc.  For  the  Voronoi 
tessellation  plot,  each  polygon  contains  a  point  and  any  point  placed  within  a  polygon  will  be 
closest  to  the  point  contained  within  that  polygon.  The  polygon  edges  are  equidistant  between 
the  points  contained  in  the  two  adjacent  polygons  and  the  triple  points  (merging  of  three  lines) 
are  equidistant  between  the  points  contained  in  the  three  adjacent  polygons.  Therefore,  the  first 
nearest  neighbors  (open  circles  in  Fig.  lc)  correspond  to  the  edges  of  the  central  polygon  (that 
contains  the  black  dot).  This  subset  of  points  is  the  maximum  number  of  nearest  neighbors  that 
the  central  point  can  have. 

In  this  manner,  the  several  techniques  can  be  modified  to  only  include  the  first  nearest  neighbors. 
For  instance,  the  modified  Wamken-Reed  method  (Fig.  Id)  only  includes  these  as  potential 
nearest  neighbors  and  cannot  expand  beyond  these,  alleviating  a  potential  problem.  Another 
method  using  these  is  to  consider  all  of  these  potential  nearest  neighbors  as  nearest  neighbors 
(Fig.  le).  Unfortunately,  this  approach  is  sensitive  to  small  perturbations  in  the  spatial  positions 
of  the  neighbors.  For  instance,  if  the  lower  right  hand  neighbor  in  Fig.  le  moves  away  from  the 
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central  point,  it  will  no  longer  share  an  edge  with  the  polygon  containing  the  black  dot;  in  this 
scenario,  the  two  adjacent  polygons  on  either  side  will  “pinch”  off  this  neighbor.  This  scenario, 
however,  has  a  physical  basis  as  these  two  dendrite  cores  will  mainly  compete  with  the  central 
core,  and  the  lower  right  core  has  a  much  less  prominent  effect  on  the  central  core.  The  last 
method,  which  will  be  examined  in  the  present  paper,  utilizes  a  criterion  based  on  the  edge 
lengths  of  the  Voronoi  polygon.  In  Figure  If,  those  neighbors  with  edge  lengths  less  than  a 
critical  fraction  of  the  total  polygon  perimeter  are  excluded  as  nearest  neighbors  (e.g.,  10%  in 
Fig.  If).  In  the  present  study,  we  consider  only  this  technique  in  evaluating  the  local  dendrite 
arm  spacing  statistics. 


<*i)  (<■)  (f) 

Figure  1.  The  difference  between  various  methods  for  defining  the  nearest  neighbors  (red  dots)  and  their  spacing 
for  a  single  point  (large  black  dot),  (a)  Initial  5x5  cubic  pattern  with  noise  fraction  of  100%  and  noise  level  of 
0.20a0.  (b)  The  Wamken-Reed  method  with  a  —  1.5  and  klmtl,i  —  3.  The  inner  circle  represents  the  average  spacing, 
dave,  of  these  neighbors  and  the  outer  circle  represents  the  cutoff  for  adding  the  next  neighbor,  dave  -hx*dsllj.  (c) 
Voronoi  tesselation  diagram  for  the  points.  The  potential  first  nearest  neighbors  are  identified  through  shared 
vertices  with  each  point,  (d)  The  modified  Wamken-Reed  method  with  a  -  1 .5  and  klmti;ii  -  3,  whereby  the  neighbors 
are  restricted  to  only  those  identified  using  the  Voronoi  tesselation  (further  explanation  of  the  difference  between  (b) 
and  (d)  is  given  in  Fig.  6).  (e)  Using  only  shared  vertices  (and  connecting  lines  forming  a  polygon)  to  identify  the 
nearest  neighbors,  (f)  Modified  tesselation-based  technique  whereby  the  nearest  neighbors  are  identified  as  those 
with  line  lengths  above  a  critical  threshold  fraction  of  the  total  perimeter  line  length  dent  of  the  tesselated  polygon 
for  the  point. 


Results 

A  micrograph  of  a  directionally-solidified  single  crystal  nickel-based  superalloy  microstructure 
that  is  polished  and  imaged  perpendicular  to  the  solidification  direction  is  shown  in  Figure  2. 
This  microstructure  was  produced  using  the  liquid  metal  cooling  technique,  as  described  in 
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Miller  (2011)  and  Elliot  et  al.  (2004).  The  black  dots  were  manually  placed  to  identify  the 
dendrite  cores.  Automated  methods  to  identify  dendrite  cores  and  other  microstructure  features 
have  been  previously  developed  and  are  invaluable  to  large  scale  analysis  (Tschopp  et  al.,  2010a, 
2010b,  2010c,  2010d).  The  white  particles  in  this  image  are  eutectic  particles.  A  total  of  393 
dendrite  cores  are  contained  in  this  image  over  an  area  of  24.25  mm2  (PDAS  of  248.4  pm).  The 
remainder  of  the  analysis  will  use  this  micrograph  as  a  template  for  characterizing  the  local 
dendrite  arm  spacing. 
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Figure  2.  Dendritic  structure  normal  to  the  withdrawal  direction  in  a  directionally-solidified  single  crystal  nickel- 
based  superalloy  cast  using  the  liquid  metal  cooling  technique  (Miller,  2010).  The  dots  denote  the  dendrite  cores. 

Edge  Detection  and  Metrics 

The  ability  to  handle  edge  effects  when  computing  local  dendrite  ann  spacings  with  dendrite 
cores  is  vital  for  quantifying  statistics  in  thin  sections,  such  as  the  wall  of  an  airfoil  blade  that 
may  only  contain  1-3  dendrite  cores  across  the  section.  As  a  first  example  of  one  such  a 
technique,  we  have  used  a  convex  hull  of  the  dendrite  cores  in  Figure  2  to  identify  “edge” 
dendrite  cores  and  quantify  the  dendrite  arm  spacing.  First,  a  convex  hull  is  generated  around 
the  points;  this  is  the  minimum  “convex”  area  that  contains  all  the  points  and  is  denoted  by  a 
blue  dotted  line  in  Figure  3.  Next,  the  edge  points  are  identified  by  finding  those  points  with 
vertices  that  lie  outside  of  the  convex  hull  (white  dots).  Then,  for  these  points,  a  new  polygon  is 
generated  by  the  intersection  of  the  initial  polygon  from  the  Voronoi  tessellation  and  the  convex 
hull  (red  polygon).  Last,  the  neighbors  can  now  be  calculated  using  either  a  new  criterion  or  the 
same  criterion  used  for  interior  points.  For  the  present  analysis,  the  same  criterion  (polygon  with 
edge  length  threshold)  was  used  for  all  points.  Further  techniques  are  needed  to  deal  with 
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complex  geometries  that  include  concave  character  and  internal  passages  in  order  to  apply  these 
techniques  to  turbine  blades. 


Figure  3.  Voronoi  tessellation  of  dendritic  structure  from  Figure  2.  The  dotted  blue  line  denotes  the  convex  hull  of 
the  dendrite  cores  and  the  red  polygons  are  those  cores  that  intersect  the  convex  hull. 

Comparison  of  Local  Metrics 

The  local  dendrite  ann  spacing  statistics  are  calculated  for  both  the  interior  and  exterior  points  to 
compare  with  the  traditional  PDAS  measurement.  The  cumulative  distribution  plot  for  the  local 
dendrite  ann  spacing  is  shown  in  Figure  4  and  the  mean  dendrite  arm  spacing  is  shown  for  both 
sites.  First,  the  interior  sites  have  a  mean  dendrite  ann  spacing  of  253.3  pm,  which  is  only  2% 
larger  than  the  bulk-measured  248.4  pm.  Moreover,  the  interior  dendrite  arm  spacings  are 
approximately  nonnally  distributed  (kurtosis  of  2.98  and  skewness  of  0.21  -  a  normal 
distribution  has  a  kurtosis  of  3  and  a  skewness  of  0).  On  the  other  hand,  the  exterior  sites  have  a 
mean  dendrite  ann  spacing  of  267.2  pm  (~5.5%  greater  than  interior  sites)  and  deviate  slightly 
from  normality  (kurtosis  of  4.49  and  skewness  of  0.47).  This  technique  allows  us  to  compare 
measurements  for  bulk  dendrite  cores  with  edge  dendrite  cores  as  well  as  measuring  the  dendrite 
ann  spacing  distribution  statistics,  which  may  be  important  for  correlating  with  both  processing 
and  properties. 

Local  Dendrite  Metrics 

The  spatial  distribution  of  local  dendrite  arm  spacing  and  coordination  number  can  provide 
insight  into  the  order/disorder  of  primary  dendrites  and  can  identify  regions  that  could  potentially 
contain  more/less  interdendritic  features  and/or  contain  different  properties.  For  instance,  the 
dendrite  ann  spacing  and  coordination  number  for  the  directionally-solidified  superalloy 
micrograph  (Fig.  2)  is  shown  in  Figure  5.  Dendrite  cores  with  local  PDAS  similar  to  the  mean 
PDAS  are  colored  white  and  those  with  PDAS  above  (below)  the  mean  PDAS  are  red  (blue).  In 
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general,  the  exterior  dendrite  cores  have  similar  PDAS  as  the  interior  dendrite  cores.  A  similar 
colorbar  is  used  for  the  coordination  number  as  well.  As  expected,  the  exterior  dendrite  cores 
have  a  lower  coordination  number  than  interior  dendrite  cores,  with  a  few  that  only  have  2 
nearest  neighbors.  However,  the  low  coordination  number  dendrite  cores  on  the  edge  are  not 
consistently  over/under  the  mean  PDAS  (i.e.,  they  do  not  significantly  bias  the  statistics  from  the 
edge  dendrite  cores). 


Figure  4.  The  cumulative  distribution  plot  of  the  primary  dendrite  arm  spacings  for  both  interior  and  exterior 
dendrite  cores  in  Figure  2.  The  solid  lines  represent  the  mean  dendrite  arm  spacings,  which  can  be  compared  to  the 
bulk  measurement  of  248.4  pm. 


Primary  dendrite  arm  spacing  (pm)  Coordination  Number 

200  220  240  260  280  300  2 _  3  4  5  _ 6 


Figure  5.  (a)  Local  dendrite  arm  spacing  (pm)  and  (b)  coordination  number  based  on  the  Voronoi  tessellation  with 
edge  length  threshold  of  10%. 
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Correlation  with  Interdendritic  Features 

The  eutectic  particles  in  Figure  2  were  segmented  to  examine  how  the  occurrence  of 
interdendritic  features  (e.g.,  pores  or  eutectic  particles)  correlates  with  local  dendrite  arm  spacing 
metrics.  The  particles  were  segmented  by  first  leveling  the  intensity  of  the  micrograph  using  a 
cubic  polynomial  with  interaction  terms  and  then  selecting  a  threshold  to  identify  the  eutectic 
particles.  Following  this  step,  the  Euclidean  distance  to  the  nearest  dendrite  core  and  the  nearest 
Voronoi  vertex  was  calculated  for  each  pixel  within  the  micrograph.  The  probability  of  a 
eutectic  particle  can  then  be  calculated  as  a  function  of  this  distance,  as  shown  in  Figure  6. 
Based  on  the  image  segmentation,  the  area  fraction  of  eutectic  particles  in  Figure  2  is  3.6%  and 
is  shown  in  red.  As  can  be  seen  from  Fig.  6a,  at  a  distance  of  approximately  100  microns,  the 
probability  of  a  eutectic  particle  exceeds  the  global  area  fraction,  i.e.,  it  is  more  favorable  for  a 
eutectic  particle  to  form.  Moreover,  this  plot  shows  that  it  is  not  favorable  for  eutectic  particles 
to  form  close  to  the  primary  dendrite  core.  Figure  6b  is  a  similar  plot  as  a  function  of  distance 
from  the  vertices  of  the  Voronoi  tessellation.  There  is  an  increased  occurrence  of  eutectic 
particles  at  vertices,  regardless  of  their  distance  from  the  dendrite  core.  While  this  analysis 
shows  the  preferential  fonnation  of  eutectic  particles  based  on  the  local  distances,  con-elation 
with  the  size  of  particles  is  also  important;  we  leave  this  for  future  work. 


(a)  (b) 

Figure  6.  The  probability  of  a  eutectic  particle  as  a  function  of  the  distance  to  (a)  the  nearest  dendrite  core  or  (b)  the 
nearest  Voronoi  vertex. 

Conclusions 

Characterizing  the  spacing  of  primary  dendrite  arms  in  directionally-solidified  microstructures  is 
an  important  step  for  developing  process-structure-property  relationships  by  enabling  the 
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quantification  of  (i)  the  influence  of  processing  on  microstructure  and  (ii)  the  influence  of 
microstructure  on  properties.  Thin-walled  directionally-solidified  structures  (e.g.,  a  turbine 
blade)  require  new  approaches  for  characterizing  the  dendrite  arm  spacing  and  the 
microstructure.  In  this  work,  we  utilized  a  new  approach  for  spatial  point  pattern  analysis  that 
was  applied  to  an  experimental  dendritic  microstructure.  This  technique  utilizes  a  Voronoi 
tessellation  of  space  surrounding  the  dendrite  cores  to  determine  nearest  neighbors  and  the  local 
dendrite  arm  spacing.  Moreover,  a  convex  hull-based  technique  was  used  to  include  edge 
effects,  which  can  be  important  for  thin  specimens.  These  methods  were  used  to  quantify  the 
distribution  of  dendrite  arm  spacings  and  their  spatial  distribution  of  dendrite  arm  spacings 
through  the  directionally-solidified  superalloy  micrograph.  Last,  the  distances  from  dendrite 
cores  and  Voronoi  vertices  were  correlated  with  the  presence  of  eutectic  particles,  which  can  be 
important  for  correlating  structure  to  mechanical  properties  for  superalloys. 
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Characterizing  the  Local  Primary  Dendrite  Arm  Spacing 
in  Directionally  Solidified  Dendritic  Microstructures 

MARK  A.  TSCHOPP,  JON  D.  MILLER,  ANDREW  L.  OPPEDAL, 
and  KIRAN  N.  SOLANKI 

Characterizing  the  spacing  of  primary  dendrite  arms  in  directionally  solidified  microstructures  is 
an  important  step  for  developing  process-structure-property  relationships  by  enabling  the 
quantification  of  (i)  the  influence  of  processing  on  microstructure  and  (ii)  the  influence  of 
microstructure  on  properties.  In  this  work,  we  utilized  a  new  Voronoi-based  approach  for 
spatial  point  pattern  analysis  that  was  applied  to  an  experimental  dendritic  microstructure.  This 
technique  utilizes  a  Voronoi  tessellation  of  space  surrounding  the  dendrite  cores  to  determine 
nearest  neighbors  and  the  local  primary  dendrite  arm  spacing.  In  addition,  we  compared  this 
technique  to  a  recent  distance-based  technique  and  a  modification  to  this  using  Voronoi  tes¬ 
sellations.  Moreover,  a  convex  hull-based  technique  was  used  to  include  edge  effects  for  such 
techniques,  which  can  be  important  for  thin  specimens.  These  methods  were  used  to  quantify  the 
distribution  of  local  primary  dendrite  arm  spacings,  their  spatial  distribution,  and  their  corre¬ 
lation  with  interdendritic  eutectic  particles  for  an  experimental  directionally  solidified  Ni-based 
superalloy  micrograph.  This  can  be  an  important  step  for  correlating  processing  and  properties 
in  directionally  solidified  dendritic  microstructures. 
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I.  INTRODUCTION 

Developing  an  enhanced  understanding  of 
mechanical  behavior  in  materials  relies  upon  sufficiently 
characterizing  microstructure  details  at  the  relevant 
length  scales  that  contribute  to  this  behavior.  Moreover, 
to  truly  enhance  the  predictive  capability  of  processing- 
structure-property  models  that  aim  to  improve  material 
performance  requires  a  quantitative  stereological 
description  of  the  relevant  microstructure  features  and, 
thereby,  the  material  itself.  Predictive  models  that 
effectively  capture  the  linkage  between  processing  and 
properties  (through  microstructure)  can  be  utilized 
within  an  integrated  computational  materials  engineer¬ 
ing  (ICME)  approach  to  design  materials  and  accelerate 
their  insertion  into  application. 

The  focus  of  the  present  work  is  on  single-crystal 
nickel-based  superalloys,  which  are  used  in  turbine 
blades  within  the  high  temperature  section  of  the 
modern  turbine  engined 1-21  In  single-crystal  nickel-based 
superalloys,  there  are  a  number  of  length  scales  of 
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microstructure  that  contribute  to  mechanical  behavior, 
ranging  from  the  /  precipitates  to  pores  and  eutectic 
particles  to  the  dendrites  themselves.  At  the  largest 
microstructure  length  scale  in  directionally  solidified 
single-crystal  microstructures,  the  features  of  interest  are 
the  dendrites;  many  features  at  lower  length  scales  (e.g., 
eutectic  particles,  precipitates,  etc.)  or  at  similar  scales 
(e.g.,  porosity,  freckle  defects,  etc.)  are  strongly  associated 
with  the  dendrite  arm  spacing  and  morphology l3-7'  The 
solidification  morphology  associated  with  dendrite  arm 
spacing  has  been  described  in  the  literature!8-9'  Histori¬ 
cally,  the  primary  dendrite  arm  spacing  (PDAS)  has  been 
found  to  correlate  with  processing  (e.g.,  solidification 
rate)17-10-14'  as  well  as  with  properties  (e.g.,  creep  strength, 
fatigue  properties)!15-6'  For  instance,  Lamm  and  Singer'6' 
produced  single-crystal  nickel-based  microstructures 
(PWA  1483)  with  a  varied  range  of  different  dendrite 
arm  spacings  (250  to  600  /mi)  and  found  that  decreasing 
the  mean  dendrite  arm  spacing  was  associated  with  an 
increased  high-cycle  fatigue  life.  The  fatigue  cracks  were 
found  to  originate  at  shrinkage  porosity  and  the  largest 
pores  correlated  with  a  large  PDAS. 

The  traditional  approach  for  measuring  PDAS  in 
single-crystal  metals,  whereby  the  number  of  dendrite 
cores  in  a  specified  area  is  related  to  the  dendrite  arm 
spacing'10-16-17'  is  given  by: 


where  2  is  PDAS,  A  is  the  area  analyzed,  n  is  the  number 
of  dendrites,  and  c  is  a  coefficient  that  depends  on  the 
microstructure.  McCartney  and  Hunt'10'  showed  that 
c  =  0.5  for  a  random  array  of  points,  c  =  1  for  a 
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square  array  of  points,  and  c  =  1.075  for  a  hexagonal 
array  of  points;  they  had  to  apply  a  correction  for  the 
bulk  dendrite  arm  spacing  2  as  processing  conditions 
caused  a  change  in  the  local  environment  of  the 
dendrites.  However,  this  approach  is  insufficient  for 
capturing  local  arm  spacings  or  the  dendrite  arm 
spacing  distribution,  and  may  provide  problems  with 
complex  geometries  such  as  turbine  blades.  In  fact,  part 
of  the  motivation  for  quantifying  the  local  PDAS  is  that 
a  narrow  distribution  (i.e.,  low  standard  deviation)  of 
local  PDAS  values  may  result  in  a  more  homogeneous 
distribution  of  interdendritic  microstructure  features 
and,  more  importantly,  a  narrow  distribution  of 
mechanical  properties.11*1  For  instance,  structural  inho¬ 
mogeneities  in  single-crystal  components,  which  may 
cause  wider  property  variations,  have  been  correlated  to 
asymmetric  heat  flux  and  transients  in  solidifications.1191 

The  research  objective  herein  is  to  evaluate  the 
capability  of  some  recent  approaches,  as  well  as  some 
modified  versions  of  these  approaches,  for  characteriz¬ 
ing  the  local  dendrite  arm  spacing  within  experimental 
dendritic  microstructures.  In  this  work,  an  experimental 
dendritic  microstructure  is  used  for  this  analysis  along 
with  three  different  techniques  that  are  based  on  the 
nearest  neighbor  spacing  and/or  a  Voronoi  tessellation 
of  the  dendrite  cores.  Comparison  of  existing  and  new 
metrics  with  traditional  PDAS  metrics  is  discussed  for 
both  local  and  global  measures.  The  current  methods 
investigated  supply  statistical  information  of  local  spac¬ 
ing  and  coordination  number  while  introducing  a 
technique  for  addressing  edge  effects  and  examining 
the  parameter  sensitivity  of  these  different  methods.  In 
comparison  to  previous  work,1201  this  work  introduces 
and  compares  the  statistical  distributions  of  local 
dendrite  arm  spacings  for  the  four  methods,  for  a  more 
quantitative  analysis.  It  was  found  that  augmenting 
existing  techniques  with  Voronoi  tessellations  to  define 
the  subset  of  first  nearest  neighbors  or  refining  existing 
Voronoi-based  techniques  resulted  in  a  more  physical 
description  of  the  local  dendrite  arm  spacing.  Moreover, 
for  certain  cases,  the  mean  local  dendrite  statistics  can 
adequately  approximate  the  PDAS  found  with  the 
traditional  bulk  characterization  technique  (Eq.  [1]). 


II.  METHODOLOGY 

The  approach  utilized  herein  to  measure  the  local 
dendrite  arm  spacing  is  based  on  a  Voronoi  tessellation 
of  the  spatial  array  of  dendrite  cores.  The  following 
analysis  techniques  were  implemented  in  MATLAB 
R2012a  (The  MathWorks.  Inc.).  To  illustrate  how  the 
present  method  works  and  differs  from  some  other 
published  methods,  we  generated  a  synthetic  5x5  cubic 
pattern  of  points  with  a  small  degree  of  noise  (100  pet 
noise  fraction,  0.20ao  noise  fraction121),  as  shown  in 
Figure  1(a).1201  For  the  purposes  of  describing  several 
different  methods  shown  in  Figure  1,  this  synthetic 
pattern  of  points  can  be  considered  as  the  cores  of 
primary  dendrites. 

One  such  method  for  measuring  the  local  dendrite 
arm  spacing  is  the  Warnken-Reed  method.121,221  The 
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Warnken-Reed  method  calculates  the  dendrite  arm 
spacing  for  a  single  point  (black  dot)  by  starting  with 
an  initial  number  of  nearest  neighbors,  A'initiai,  (3  closest 
neighbors)  and  iteratively  adding  potential  nearest 
neighbors  that  are  within  a  cutoff  distance  defined  by 
the  already-added  nearest  neighbors.  For  instance,  the 
inner  circle  in  Figure  1(b)  represents  the  average  spac¬ 
ing,  rfavg,  of  these  neighbors  and  the  outer  circle  repre¬ 
sents  the  cutoff  for  adding  the  next  neighbor,  r/avg  + 
xrfsld,  where  rfsld  is  the  standard  deviation  of  the  nearest 
neighbor  spacings  and  a  is  a  parameter  that  is  typically 
between  1  and  2.  Neighbors  continue  to  be  added  until 
the  cutoff  does  not  include  any  new  neighbors.  The  local 
coordination  number  and  dendrite  arm  spacing  is 
calculated  from  the  neighbors  added  (shown  as  red 
dots).  However,  if  the  standard  deviation  of  the  dis¬ 
tances  of  the  nearest  neighbors  r/sld  or  the  parameter  x  is 
large,  this  technique  can  continue  to  add  nearest 
neighbors  beyond  the  first  nearest  neighbors;  our 
implementation  stopped  after  20  nearest  neighbors. 
Clearly,  a  method  for  restricting  the  number  of  nearest 
neighbors  using  such  a  technique  is  necessary. 

A  simple  way  of  identifying  the  potential  first  nearest 
neighbors  is  to  perform  a  Voronoi  tessellation  of  the 
space  surrounding  the  points,  as  shown  in  Figure  1(c). 
The  polygon  edges  are  equidistant  between  the  points 
contained  in  the  two  adjacent  polygons  and  the  triple 
points  (merging  of  three  lines)  are  equidistant  between 
the  points  contained  in  the  three  adjacent  polygons. 
Therefore,  the  first  nearest  neighbors  (FNNs,  shown  as 
open  circles  in  Figure  1(c))  correspond  to  the  edges  of 
the  central  polygon  (that  contains  the  black  dot).  This 
subset  of  points  is  the  maximum  number  of  nearest 
neighbors  that  the  central  point  can  have. 

In  this  manner,  several  techniques  have  been  identi¬ 
fied  to  quantify  a  local  dendrite  arm  spacing  based  on 
the  Voronoi-identified  FNNs.1201  For  instance,  the 
Voronoi  Warnken-Reed  method  (Figure  1(d))  only 
includes  the  Voronoi  FNNs  as  potential  nearest  neigh¬ 
bors  and  cannot  expand  beyond  these,  alleviating  a 
potential  problem  of  selecting  second  nearest  neighbors 
or  greater.  Another  method  using  the  Voronoi  FNNs  is 
to  consider  all  of  these  potential  nearest  neighbors  as 
nearest  neighbors  (Figure  1(e)),  as  in  Brundidge  el  a/.1141 
Unfortunately,  this  approach  is  sensitive  to  small 
perturbations  in  the  spatial  positions  of  the  neighbors. 
For  instance,  if  the  lower  right  hand  neighbor  in 
Figure  1(e)  moves  away  from  the  central  point,  it  no 
longer  shares  an  edge  with  the  polygon  containing  the 
black  dot;  in  this  scenario,  the  two  adjacent  polygons  on 
either  side  effectively  “pinch”  off  this  neighbor.  This 
scenario,  however,  has  a  physical  basis  as  these  two 
dendrite  cores  mainly  compete  with  the  central  core,  and 
the  lower  right  core  has  a  much  less  prominent  effect  on 
the  central  core.  The  last  method,  which  is  examined  in 
the  present  paper,  utilizes  a  criterion  based  on  the  edge 
lengths  of  the  Voronoi  polygon.  In  Figure  1(f),  those 
neighbors  with  edge  lengths  less  than  a  critical  fraction, 
i/crit,  of  the  total  polygon  perimeter  are  excluded  as 
nearest  neighbors  ( e.g .,  10  pet  in  Figure  1(f)).  In 
the  present  study,  the  local  dendrite  arm  spacing 
statistics  are  evaluated  using  these  four  techniques: 
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Fig.  1 — The  difference  between  various  methods  for  defining  the  nearest  neighbors  (red  dots)  and  their  spacing  for  a  single  point  (large  black 
dot),  (a)  Initial  5x5  cubic  pattern  with  noise  fraction  of  100  pet  and  noise  level  of  0.20uo-  (b)  The  Warnken-Reed  method  with  tz  =  1.5  and 
^initial  =  3-  The  inner  circle  represents  the  average  spacing,  davs,  of  these  neighbors  and  the  outer  circle  represents  the  cutoff  for  adding  the  next 
neighbor,  r/avg  +  v£/sti1.  (c)  Voronoi  tessellation  diagram  for  the  points.  The  potential  first  nearest  neighbors  are  identified  through  shared  vertices 
with  each  point,  (d)  The  modified  Warnken-Reed  method  with  a  =  1.5  and  kinitia j  =  3,  whereby  the  neighbors  are  restricted  to  only  those  iden¬ 
tified  using  the  Voronoi  tessellation,  (e)  Using  only  shared  vertices  (and  connecting  lines  forming  a  polygon)  of  the  Voronoi  tessellation  to  iden¬ 
tify  the  nearest  neighbors  (dcrit  =  0.0).  (/)  Modified  tessellation-based  technique  whereby  the  nearest  neighbors  are  identified  as  those  with  line 
lengths  above  a  critical  threshold  fraction  of  the  total  perimeter  line  length  t/cr jt  =  0.10  of  the  tessellated  polygon  for  the  point  (Reprinted 
from^)  (Color  figure  online). 


Warnken-Reed,  Voronoi  Warnken-Reed,  and  the 
Voronoi  technique  with  (c?crit  >  0)  and  without 
(^erit  =  0)  a  line  length  threshold. 

As  an  example  of  a  more  disordered  structure, 
Figure  2  plots  the  four  different  methods  for  a  different 
configuration  of  surrounding  points  (dendrite  cores).  In 
Figure  2(b),  the  iterative  Warnken-Reed  method  con¬ 
tinues  to  non-physically  add  neighbors  beyond  the  first 
nearest  neighbors  due  to  a  large  initial  k/std  value  from 
the  initial  three  distances.  The  Voronoi-modified  version 
in  Figure  2(d)  stops  at  four  nearest  neighbors  despite  the 
fact  that  several  points  lie  within  the  outer  boundary 
computed  by  this  method.  The  Voronoi  method  with 
dent  =  0.0  clearly  overestimates  the  number  of  nearest 
neighbors,  while  the  four  nearest  neighbors  identified 
through  dCIit  =  0.10  (Figure  2(f))  perhaps  offers  a  better 
approximation  of  the  number  of  nearest  neighbors. 
Interestingly,  comparing  the  methods  in  Figures  2(d) 
and  (f),  the  coordination  number  is  the  same,  but  the 
nearest  neighbors  identified  is  different.  This  is  due  to 
the  Warnken-Reed  method  being  a  distance-based 
method,  and  identifying  the  four  closest  neighbors, 
while  the  modified  Voronoi  technique  is  based  on  the 
edge  lengths  of  the  Voronoi  polygon,  and  hence  utilizes 


this  to  identify  nearest  neighbors  (which  may  not  be  the 
closest  neighbors). 

The  traditional  PDAS  metric  does  not  consider  the 
order  or  disorder  of  the  dendrites  within  the  micro¬ 
structure.  Figure  2  illustrates  why  a  local  metric  for 
PDAS  may  be  needed.  For  the  field  of  view  given  in 
Figures  1(a)  and  2(a),  the  bulk  PDAS  metric  would  be 
the  same  since  the  number  of  dendrites  n  and  the  area  A 
are  equal  (see  Eq.  [1]).  However,  the  disorder  of  the 
dendritic  structure  in  the  case  of  Figure  2  may  yield  (i)  a 
more  uneven  distribution  of  solute  elements,  (ii)  the 
formation  of  second  phase  particles,  (iii)  the  formation 
of  gas  or  shrinkage  porosity,  or  (iv)  the  lateral  growth  of 
secondary  dendrite  arms.  Hence,  in  addition  to  the  bulk 
PDAS  values,  understanding  how  processing  conditions 
may  impact  the  disorder  of  the  dendritic  structure  may 
be  important  for  understanding  the  properties  of  direc¬ 
tionally  solidified  alloys. 

Other  techniques  exist  for  quantifying  the  homoge¬ 
neity  or  heterogeneity  of  PDAS  in  directionally  solidi¬ 
fied  dendritic  microstructures.  For  instance,  the  minimal 
spanning  tree  (MST)  method12'1  provides  a  statistical 
analysis  of  the  disorder  in  a  system  of  points  by 
connecting  all  points  with  the  shortest  line  segments 
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Fig.  2 — The  difference  between  various  methods  for  defining  the  nearest  neighbors  (red  dots)  and  their  spacing  for  a  single  point  (large  black 
dot)  with  a  distorted  local  environment.  Parts  ( a )  through  (/)  are  as  Fig.  1:  (a)  initial  pattern,  (b)  the  Warnken-Reed  method  with  a  =  1.5  and 
initial  =  3,  (c)  the  Voronoi  tessellation  diagram  for  the  points,  (d)  the  Voronoi-modified  Warnken-Reed  method  with  a  =  1.5  and  &initia i  =  3, 
(e)  the  Voronoi  method  (dct it  =  0.0),  and  (f)  the  modified  tessellation-based  technique  with  dent  =  0.10  (Color  figure  online). 


possible.  In  this  manner,  the  mean  distance  of  all  line 
segments  ( m )  and  the  standard  deviation  (rr)  characterize 
the  disorder  of  the  system  and  casting  these  values  into  a 
m-tr  design  space  allows  for  comparison  between  differ¬ 
ent  point  systems.12  '1  This  has  been  effectively  applied  to 
characterize  the  mean  dendrite  arm  spacing,  PDAS 
distribution,  and  the  disorder  in  first  Pb-Tl  alloys1241  and 
subsequently  in  other  alloy  systems. ll1,18’251  As  an 
example  of  this  technique.  Figure  3  plots  the  dendrite 
cores  and  connecting  line  segments  for  the  single-crystal 
nickel-based  superalloy  micrograph  used  in  this  study 
(Figure  4).  Moreover,  other  methods  such  as  radial 
distribution  functions,  fast  Fourier  transforms,  and/or 
correlation  functions  can  also  be  used  to  characterize  the 
dendrite  arm  spacing  distribution.  However,  it  should  be 
noted  that  these  approaches  are  not  intended  for  local 
characterization  of  the  dendrite  arm  spacing  and  are  not 
as  effective  for  correlating  the  local  spacing  with  local 
microstructure  features  as  shown  herein.  Moreover, 
these  techniques  do  not  quantify  the  number  of  nearest 
neighbors  and  are  often  coupled  with  Voronoi  polygons 
to  compute  the  nearest  neighbor  distributions.  Rather, 
these  analysis  methods  are  more  effective  at  character¬ 
izing  and  comparing  the  homogeneity/heterogeneity  of 
the  dendritic  structure  between  different  processing 
conditions.  Hence,  there  will  be  limited  discussion  of 
these  techniques  in  the  present  work. 
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III.  RESULTS 

A.  Application  to  Dendritic  Microstructure 

A  micrograph  of  a  directionally  solidified  single¬ 
crystal  nickel-based  superalloy  microstructure  that  is 
polished  and  imaged  perpendicular  to  the  solidification 
direction  is  shown  in  Figure  4.  This  microstructure  was 
produced  using  the  liquid  metal  cooling  technique,  as 
described  in  Miller1261  and  Elliott  et  all4 11  First,  the 
dendrite  cores  were  identified  manually  (white  and  black 
dots).  Automated  methods  to  identify  dendrite  cores  can 
be  invaluable  for  future  large  scale  analysis.127’281  More¬ 
over,  the  white  particles  in  this  image  are  eutectic 
particles.  A  total  of  393  dendrite  cores  are  contained  in 
this  image  over  an  area  of  24.25  mm2,  giving  a  PDAS  of 
248.4  pm  using  c  =  1  (Eq.  [1]).  The  remainder  of  the 
analysis  uses  this  micrograph  as  a  template  for  charac¬ 
terizing  the  local  dendrite  arm  spacing. 

B.  Accounting  for  Image) Part  Edge  Effects 

The  ability  to  handle  edge  effects  when  computing 
local  dendrite  arm  spacings  with  dendrite  cores  is  vital 
for  quantifying  statistics  in  thin  sections,  such  as  the 
wall  of  an  airfoil  blade  that  may  only  contain  1  to  3 
dendrite  cores  across  the  section/"7’281  As  a  first  example 
of  one  such  a  technique,  we  have  used  a  convex  hull  of 
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Fig.  3 — Minimal  spanning  tree  method^23'  for  defining  the  spacing 
and  homogeneity  of  a  microstructure  (set  of  points),  whereby  the 
lines  represent  the  minimal  distance  of  connecting  line  segments.  The 
set  of  points  selected  for  this  example  were  selected  from  the  den¬ 
drite  cores  shown  in  Fig.  4,  where  the  white  dots  indicate  ‘edge’  den¬ 
drite  cores. 


Fig.  4 — Dendritic  structure  normal  to  the  withdrawal  direction  in  a 
directionally  solidified  single-crystal  nickel-based  superalloy  cast 
using  the  liquid  metal  cooling  technique)26^  The  dots  denote  the  den¬ 
drite  cores,  where  the  white  dots  indicate  ‘edge’  dendrite  cores,  as 
discussed  in  Fig.  5  and  the  associated  text. 

the  dendrite  cores  in  Figure  4  to  identify  “edge” 
dendrite  cores  and  quantify  the  dendrite  arm  spacing. 
The  dendrite  core  locations  are  first  extracted  from  the 
experimental  image,  as  shown  in  Figure  5.  Then,  a 
convex  hull  is  generated  around  the  points;  this  is  the 
minimum  "convex”  area  that  contains  all  the  points. 
Next,  the  edge  points  (white  dots  in  Figure  4)  are 
identified  by  finding  those  points  with  Voronoi  vertices 
that  lie  outside  of  the  convex  hull  (dotted  blue  line  in 
Figure  5(a)).  Then,  to  utilize  Voronoi-based  techniques 
for  these  points,  a  new  polygon  is  generated  by  the 
intersection  of  the  initial  polygon  from  the  Voronoi 
tessellation  and  the  convex  hull;  the  new  polygon  of  the 
edge  dendrite  cores  is  colored  red  in  Figure  5(a)  to 
distinguish  from  the  bulk  dendrite  cores.  The  polygons 
belonging  to  the  interior  and  edge  dendrites  are  shown 
in  Figures  5(b)  and  (c),  with  a  random  coloring  scheme 


(a)  Convex  hull 


(b)  Interior  dendrite  cores  (c)  Edge  dendrite  cores 

Fig.  5 — (a)  Voronoi  tessellation  of  dendritic  structure  from  Fig.  4. 
The  dotted  blue  line  (surrounding  the  points)  denotes  the  convex 
hull  of  the  dendrite  cores  and  the  red  polygons  delineate  the  cores 
that  intersect  the  convex  hull.  The  interior  and  edge  dendrites  are 
shown  in  ( b )  and  (c),  respectively,  with  each  polygon  colored  differ¬ 
ently  as  a  guide  to  the  eye  (Color  figure  online). 

used  to  delineate  the  different  polygons.  Last,  the 
neighbors  can  now  be  calculated  using  either  a  new 
criterion  or  the  same  criterion  used  for  interior  points. 
For  the  present  analysis,  the  same  criterion  (polygon 
with  edge  length  threshold)  was  used  for  all  points; 
although  herein  the  interior  dendrite  cores  are  used  to 
compare  statistics  with  other  techniques  and  bulk  PDAS 
values.  More  complicated  techniques  are  needed  to  deal 
with  complex  geometries  that  include  concave  character 
and  internal  passages  in  order  to  eventually  apply  these 
techniques  to  complex  structures  such  as  turbine  blades. 
Multiple  instantiations  of  microstructures  with  edge 
effects  can  shed  light  on  the  appropriate  method  for 
determining  the  local  PDAS  at  edges,  which  may  be 
different  from  that  used  in  the  interior. 

C.  Spatial  Distribution  of  Local  Primary  Dendrite 
Arm  Spacings 

The  spatial  distribution  of  local  dendrite  arm  spacing 
and  coordination  number  can  provide  insight  into  the 
order/disorder  of  primary  dendrites  and  can  identify 
regions  that  could  potentially  contain  more/less  inter- 
dendritic  features  and/or  contain  different  properties. 
For  instance,  the  PDAS  and  coordination  number  for  the 
directionally  solidified  superalloy  micrograph  (Figure  4) 
is  shown  in  Figure  6.  In  this  example,  we  used  the 
Voronoi  tessellation-based  technique  with  an  edge  length 
threshold  of  dciit  =  0.12.  Dendrite  cores  with  local 


430— VOLUME  45A,  JANUARY  2014  METALLURGICAL  AND  MATERIALS  TRANSACTIONS  A 


28 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


PDAS  similar  to  the  mean  PDAS  of  the  bulk  (248.4  /tm) 
are  colored  white  and  those  with  PDAS  above  (below) 
the  mean  PDAS  are  red  (blue);  the  lower  and  upper 
bounds  of  the  color  bar  are  —25  and  +  25  pet  of  the  mean 
PDAS  value,  respectively.  In  general,  the  exterior  den- 


200  220  240  260  280  300 


(a)  Primary  dendrite  arm  spacing  (fim) 
2  3  4  5  6 


Fig.  6 — ( a )  Local  dendrite  arm  spacing  (/zm)  and  ( b )  coordination 
number  based  on  the  Voronoi  tessellation  with  edge  length  threshold 
of  c/crit  =  0.12  or  12  pet. 


drite  cores  have  similar  PDAS  as  the  interior  dendrite 
cores  using  this  technique.  A  similar  color  bar  is  used  for 
the  coordination  number  as  well.  As  would  be  expected, 
the  exterior  dendrite  cores  tend  to  have  a  lower  coordi¬ 
nation  number  than  the  interior  dendrite  cores,  with  a 
few  that  only  have  2  nearest  neighbors.  However,  the 
dendrite  cores  with  a  low  coordination  number  on  the 
edges  are  not  consistently  over/under  the  mean  PDAS 
(;>.,  they  do  not  significantly  bias  the  statistics  from  the 
edge  dendrite  cores).  Future  work  will  examine  what 
techniques  may  be  most  applicable  for  characterizing 
local  dendrite  arm  spacings  and  coordination  numbers 
for  dendrite  cores  on  free  surfaces.  It  is  envisioned  that 
sectioning  large  numbers  of  instantiations  of  syntheti¬ 
cally  generated  microstructures  of  known  bulk  dendrite 
arm  spacings  can  be  used  to  understand  the  bias 
introduced  by  edge  effects  and  to  understand  what  are 
the  best  techniques  for  quantifying  the  local  spacing. 

The  local  dendrite  arm  spacing  for  the  remaining  three 
techniques  is  shown  in  Figure  7.  The  same  color  bar  for 
local  PDAS  as  in  Figure  6  is  used  here.  First,  notice  that 
the  Voronoi  tessellation-based  technique  with  an  edge 
length  threshold  of  dc rit  =  0.0  has  a  much  larger  fraction 
of  dendrite  cores  with  PDAS  greater  than  the  bulk  PDAS 
than  below  the  bulk  PDAS  (83.5  pet  above  248.4  /mi). 
Clearly,  the  local  PDAS  is  overpredicted  in  this  case.  The 
Warnken-Reed  and  Voronoi  Warnken-Reed  methods 
are  shown  in  Figures  7(b)  and  (c).  At  first  glance,  a 
majority  of  the  local  PDAS  values  are  very  similar  between 
the  two  methods  (-79  pet).  However,  -21  pet  of  the  cores 
resulted  in  a  difference  between  the  two  techniques,  which 
is  caused  by  the  original  Warnken-Reed  method  using 
neighbors  outside  of  those  FNNs  identified  from  the 
Voronoi  polygons.  In  every  case,  the  Warnken-Reed 
method  resulted  in  higher  local  PDAS  values  than  the 
Voronoi  Warnken-Reed  method,  as  would  be  expected 
since  this  is  purely  a  distance-based  method  and  sub¬ 
sequent  additions  can  only  increase  the  local  PDAS. 

Figure  8  shows  the  difference  in  local  PDAS  values 
between  the  two  techniques.  In  several  cases,  the 
difference  is  greater  than  250  /tm  and/or  100  pet  of  the 
PDAS  value  quantified  by  the  Voronoi  Warnken-Reed 
method.  The  differing  dendrite  cores  is  approximately 


200  220  240  260  280  300 


(a)  Voronoi  Tessellation  ( dcri ,  =  0.0)  (b)  Warnken-Reed  (a  =  2.0)  (c)  Voronoi  Warnken-Reed  (a  =  2.0) 

Fig.  7 — Local  dendrite  arm  spacing  (^m)  for  the  three  techniques  not  shown  in  Fig.  6:  ( a )  Voronoi  tessellation  with  edge  length  threshold  of 
dCTit  =  0.0,  ( b )  Warnken-Reed  technique  with  a  =  2.0,  and  (c)  Voronoi  Warnken-Reed  with  a  =  2.0. 
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an  equal  percentage  for  edge  dendrites  as  well  as  interior 
dendrite  cores. 

D.  Local  Primary  Dendrite  Arm  Spacing  Statistics 

The  local  dendrite  arm  spacing  statistics  are  also 
calculated  for  the  interior  dendrite  cores  to  compare 
with  the  traditional  PDAS  measurement.  The  cumula¬ 
tive  distribution  function  plot  for  the  local  dendrite  arm 
spacing  is  shown  in  Figure  9  for  the  three  different 
techniques  over  a  range  of  parameter  values,  which  are 


0  50  100  150  200  250 


Primary  dendrite  arm  spacing  ( pm)  difference 


Fig.  8 — Difference  in  the  local  dendrite  arm  spacing  (pm)  between 
the  Warnken-Reed  and  Voronoi  Warnken-Reed  techniques  with 
a  =  2.0.  The  Warnken-Reed  technique  had  a  greater  PDAS  value  in 
every  case  (—21  pet  of  dendrite  cores  are  different). 


given  in  the  legend.  The  bulk  PDAS  measurement  is 
shown  as  a  vertical  black  line  and  the  hexagonal  star 
shows  the  50th  percentile  intersection  point.  The  local 
dendrite  arm  spacing  distributions  are  characterized  in 
terms  of  mean,  standard  deviation,  skewness,  and 
kurtosis  (Table  I),  while  the  coordination  number  dis¬ 
tributions  are  characterized  in  terms  of  their  mean  and 
the  percentages  of  3,  4,  5,  6,  and  7  +  nearest  neighbors 
(Table  II).  The  skewness  and  kurtosis  measure  the 
asymmetry  of  the  distribution  and  the  sharpness  of  the 
peak/thickness  of  the  tail,  respectively.  The  skewness 
and  kurtosis  are  0  and  3,  respectively,  for  a  normal 
distribution. 

There  are  distinct  differences  between  the  local 
dendrite  arm  spacing  distributions  calculated  by  the 
four  techniques  (Figure  9  and  Table  I).  The  Warnken- 
Reed  and  Voronoi  Warnken-Reed  are  compared  ini¬ 
tially.  In  the  case  of  the  Warnken-Reed  method,  the 
PDAS  distribution  is  shifted  towards  large  PDAS  values 
at  high  tz  values  (a  positive  skewness  value  gives  a  long 
tail)  and  has  a  sharper  peak  and  a  longer,  fatter  tail 
(high  kurtosis  values),  more  so  than  the  other  methods. 
This  skewness  is  caused  by  an  overestimation  of  the 
number  of  nearest  neighbors  in  some  cases,  due  to  large 
values  of  either  dsid  or  the  parameter  a.  For  the  Voronoi 
Warnken-Reed  method,  there  is  a  lack  of  a  long  tail  and 
the  skewness/kurtosis  of  the  distribution  tends  more 
towards  normality.  However,  the  calculated  mean 
PDAS  with  this  method  tends  to  underestimate  the 
bulk-measured  PDAS.  While  the  maximum  number  of 
nearest  neighbors  (8  for  a  >  1.2)  is  more  realistic,  a  large 
percentage  of  dendrite  cores  are  predicted  to  have  only  3 
nearest  neighbors,  even  in  the  case  of  a  large  a 
parameter  (48.6  pet  for  a  =  2.0).  It  is  also  interesting 
that  increasing  the  a  parameter  for  the  case  of  the 


Local  Primary  Dendrite  Arm  Spacing  (  pm)  Local  Primary  Dendrite  Arm  Spacing  (  pm) 


(a) 


(b) 


Fig.  9 — Probability  distribution  functions  for  the  various  local  characterization  methods  compared  within  for  the  internal  dendrites  within  the 
dendritic  microstructure  shown  in  Fig.  4.  The  four  different  techniques  are  compared  with  the  bulk  PDAS  for  a  range  of  parameter  values.  The 
upper  bound  of  the  parameter  range  for  each  technique  is  shown  as  a  dotted  line.  To  facilitate  the  comparison,  the  Warnken-Reed  and  the 
Voronoi  technique  (dcrit  =  0.0)  are  shown  in  (a),  and  the  remaining  Voronoi -modified  techniques  are  shown  in  ( b ). 
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Table  I.  Local  Primary  Dendrite  Arm  Spacing  Statistics 


Method 

Parameter 

Primary  Dendrite  Arm  Spacing 

Mean  (/mi) 

Diff.  (Pet) 

Std  (jum) 

Skewness 

Kurtosis 

Bulk  (Eq.  [1],  c  =  1) 

— 

248.4 

0.0 

— 

— 

— 

Voronoi  Tessellation 

dcrit  =  0.00 

272.9 

9.9 

28.0 

0.1 

3.3 

Voronoi  Tessellation 

dcrit  =  0.02 

270.0 

8.7 

27.2 

0.1 

3.2 

Voronoi  Tessellation 

dcrit  ~  0-04 

266.4 

7.2 

26.8 

0.1 

3.1 

Voronoi  Tessellation 

dcrit  =  0.06 

263.0 

5.9 

25.7 

0.1 

3.1 

Voronoi  Tessellation 

deni  =  0.08 

258.4 

4.0 

26.1 

0.2 

3.1 

Voronoi  Tessellation 

(4rit  =  0.10 

253.3 

2.0 

25.3 

0.2 

3.0 

Voronoi  Tessellation 

dcrit  =  0.12 

247.6 

-0.3 

26.0 

0.2 

2.9 

Voronoi  Warnken-Reed 

a  =  1.0 

230.0 

-7.4 

25.0 

0.4 

3.5 

Voronoi  Warnken-Reed 

a  =  1.2 

230.9 

-7.0 

26.0 

0.4 

3.4 

Voronoi  Warnken-Reed 

a  =  1.4 

232.7 

-6.3 

27.4 

0.4 

3.2 

Voronoi  Warnken-Reed 

a  =  1.6 

236.1 

-5.0 

29.5 

0.4 

3.1 

Voronoi  Warnken-Reed 

a  =  1.8 

239.0 

-3.8 

30.5 

0.4 

3.0 

Voronoi  Warnken-Reed 

a  =  2.0 

242.2 

-2.5 

31.8 

0.3 

3.0 

W  arnken-Reed 

a  =  1.0 

230.1 

-7.4 

25.1 

0.3 

3.4 

W  arnken-Reed 

a  =  1.2 

231.8 

-6.7 

26.9 

0.5 

3.5 

W  arnken-Reed 

a  =  1.4 

234.5 

-5.6 

29.4 

0.6 

3.8 

W  arnken-Reed 

a  -  1.6 

239.2 

-3.7 

33.0 

0.8 

4.0 

W  arnken-Reed 

a  =  1.8 

248.1 

-0.1 

48.1 

1.9 

8.3 

W  arnken-Reed 

a  =  2.0 

259.2 

4.3 

64.2 

1.9 

6.5 

Minimal  Spanning  Tree 

N/A 

215.2 

-13.4 

34.1 

-0.5 

4.8 

Table  II. 

Local  Coordination  Number  Statistics 

Method 

Parameter 

3 

Coordination  Number  (Pet) 

4  5  6 

>7 

Mean 

Voronoi  Tessellation 

dcrit  =  0.00 

0.0 

2.5 

20.4 

57.3 

19.8 

5.98 

Voronoi  Tessellation 

dcrit  =  0-02 

0.0 

4.0 

26.0 

57.0 

13.0 

5.80 

Voronoi  Tessellation 

dcrit  =  0.04 

0.0 

7.4 

34.7 

51.4 

6.5 

5.57 

Voronoi  Tessellation 

dcrit  =  0.06 

0.0 

11.8 

45.2 

39.9 

3.1 

5.35 

Voronoi  Tessellation 

dcrit  =  0.08 

1.5 

24.1 

49.2 

24.8 

0.3 

4.98 

Voronoi  Tessellation 

dcrit  =  0.10 

5.3 

43.0 

41.5 

10.2 

0.0 

4.57 

Voronoi  Tessellation 

dcrit  =  0.12 

17.3 

52.3 

26.9 

3.4 

0.0 

4.16 

Voronoi  Warnken-Reed 

or  =  1.0 

94.1 

4.6 

0.9 

0.0 

0.3 

3.08 

Voronoi  Warnken-Reed 

a  =  1.2 

87.6 

8.7 

2.8 

0.3 

0.6 

3.18 

Voronoi  Warnken-Reed 

a  =  1.4 

79.3 

11.8 

5.3 

2.5 

1.2 

3.35 

Voronoi  Warnken-Reed 

a  =  1.6 

65.9 

18.0 

8.0 

5.3 

2.8 

3.62 

Voronoi  Warnken-Reed 

or  =  1.8 

56.3 

19.8 

12.1 

8.4 

3.4 

3.84 

Voronoi  Warnken-Reed 

a  =  2.0 

48.6 

19.8 

14.9 

11.5 

5.3 

4.07 

W  arnken-Reed 

a  =  1.0 

91.6 

7.1 

0.9 

0.0 

0.3 

3.10 

Warnken-Reed 

or  =  1.2 

83.9 

9.6 

4.0 

1.2 

1.2 

3.27 

Warnken-Reed 

a  =  1.4 

72.8 

14.6 

7.1 

2.8 

2.8 

3.52 

W  arnken-Reed 

a  —  1.6 

58.5 

21.7 

9.0 

4.0 

6.8 

3.89 

W  arnken-Reed 

a  =  1.8 

52.0 

19.2 

10.8 

6.2 

11.8 

4.63 

W  arnken-Reed 

a  =  2.0 

44.0 

19.2 

12.1 

7.1 

17.6 

5.55 

Voronoi  Warnken-Reed  method  tends  to  shift  the  slope 
of  the  probability  distribution  function  without  affecting 
either  the  minimum  or  maximum  local  dendrite  arm 
spacings. 

For  comparison,  the  minimal  spanning  tree  method 
(Figure  3)  was  also  included  in  Table  I.  Not  surpris¬ 
ingly,  the  mean  distance  of  the  connecting  line  segments 
is  much  shorter  than  the  bulk  calculated  PDAS  using 
Eq.  [1]  with  c  —  1.  Remember  that  the  MST  method  is 
composed  of  the  shortest  line  segments  to  connect  all 
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dendrites.  Both  the  MST  standard  deviation  and  kur- 
tosis  values  are  larger  than  the  Voronoi  tessellation 
method  (for  all  dcr *)  and  the  Voronoi  Warnken-Reed 
method  (for  all  a),  indicating  a  wider  distribution  and  a 
larger  deviation  of  the  distribution  from  normality 
(kurtosis  =  3).  Moreover,  the  distribution  is  skewed 
towards  a  larger  tail  at  lower  distances  (negative 
skewness)  unlike  the  other  techniques,  which  again  is 
associated  with  the  selection  of  the  shortest  line  seg¬ 
ments  to  characterize  the  spacing. 
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Quantifying  the  coordination  number  and  the  local 
PDAS  values  using  all  FNNs  identified  by  the  Voronoi 
tessellation  polygons  (dCI jt  =  0)  clearly  overestimates 
both  measures;  mean  PDAS  is  -10  pet  off  from  the 
bulk-measured  PDAS  value  and  -20  pet  of  dendrite 
cores  have  more  than  6  nearest  neighbors.  As  the  edge 
length  threshold  parameter  increases,  less  nearest  neigh¬ 
bors  are  identified  and  the  calculated  mean  PDAS 
approaches  the  bulk-measured  PDAS  value  (within  0.3 
pet  for  dail  =  0.12).  For  t/orit  =  0.12,  the  majority  of 
dendrite  cores  have  4  nearest  neighbors  (>50  pet), 
followed  by  5  nearest  neighbors  (26.9  pet)  and  3  nearest 
neighbors  (17.3  pet).  Moreover,  the  local  PDAS  distri¬ 
bution  has  a  low  skewness  value  (0.2)  and  a  kurtosis  of 
2.9,  indicating  an  approximately  normal  distribution.  In 
general,  the  Voronoi  tessellation-based  technique  with 
an  edge  length  threshold  criterion  of  dCI jt  =  0.12 
tends  to  give  the  best  agreement  in  terms  of  both 
bulk-measured  PDAS  and  coordination  number. 

E.  Correlation  with  Interdendritic  Features 

The  relationship  between  the  occurrence  of  interden¬ 
dritic  features  ( e.g .,  pores  or  eutectic  particles)  and  the 
local  dendrite  arm  spacing  (or  distance  from  cores,  etc.) 
can  provide  insight  into  the  importance  of  quantifying 
the  local  spacings.  We  have  examined  how  these 
metrics  may  relate  to  the  formation  of  eutectic  particles 
in  this  work  by  first  segmenting  the  interdendritic 
particles  and  then  computing  probability  distribution 
functions. 

The  eutectic  particles  in  Figure  4  were  segmented 
using  the  following  process.  The  particles  were  seg¬ 
mented  by  first  leveling  the  intensity  of  the  micrograph 
using  a  cubic  polynomial  with  interaction  terms.  This 
step  ensured  that  there  was  no  shift  in  contrast  from  one 
side  of  the  micrograph  to  the  other  (due  to  uneven 
etching,  etc.).  The  threshold  intensity  was  then  selected 
by  maximizing  the  difference  in  the  mean  intensity 
between  the  two  distributions  (eutectic  particle 
and  matrix).  Then,  a  size  threshold  was  enforced  by 


discarding  eutectic  particles  with  less  than  5  pixels  (1  pixel 
-1.7  /tm,  i.e.,  5  pixels  =  15.2  /an2).  As  an  example  of  the 
segmentation.  Figure  10(a)  shows  a  1mm  x  1  mm  region 
from  Figures  4  and  10(b)  shows  the  corresponding  binary 
image  of  the  segmented  eutectic  particles  (in  white).  While 
other  segmentation  techniques,  such  as  region  grow¬ 
ing,12^  have  been  used  to  segment  these  interdendritic 
particles  features  in  3D  serial  slices,  the  method  utilized 
within  produced  sufficient  segmentation  of  the  particles 
for  the  subsequent  analysis. 

The  Euclidean  distance  to  the  nearest  dendrite  core 
and  the  nearest  Voronoi  vertex  was  then  calculated  for 
each  pixel  within  the  micrograph.  The  Euclidean  dis¬ 
tance  is  the  distance  from  each  pixel  to  the  nearest 
feature,  which  in  this  case  is  either  the  centroids  of  the 
dendrite  cores  or  the  Voronoi  vertices,  and  this  metric  is 
repeated  over  all  pixels  within  the  image  to  create  a  map. 
As  an  example.  Figure  10(c)  shows  the  Euclidean 
distance  map  for  the  same  1  mm  x  1  mm  area  utilizing 
the  dendrite  core  centroids  identified  in  Figure  4.  The 
darker  intensity  indicates  closer  Euclidean  distances  to 
the  dendrite  core  and  the  lightest  pixels  between  the 
dendrite  cores  actually  correspond  to  the  boundaries  of 
the  Voronoi  tessellation. 

The  probability  of  encountering  (or  forming)  a  eutectic 
particle  can  then  be  calculated  as  a  function  of  this 
Euclidean  distance  from  the  nearest  dendrite  core  or  the 
Voronoi  vertex,  as  shown  in  Figure  11.  Based  on  the 
image  segmentation,  the  area  fraction  of  eutectic  particles 
in  Figure  4  is  3.6  pet  and  is  shown  as  a  red  line  in 
Figure  1 1.  The  pixels  lying  within  100  /tm  of  the  image 
boundaries  were  excluded  to  eliminate  the  possibility  that 
dendrite  cores  just  outside  of  the  image  field  of  view  could 
affect  the  statistics.  As  can  be  seen  from  Figure  11(a),  the 
left  (right)  blue  line  indicates  the  distance  whereby  all 
distances  below  (above)  have  a  probability  of  having  a 
eutectic  particle  that  is  lower  (higher)  than  the  global  area 
fraction  (red  line),  i.e..  it  is  less  (more)  favorable  for  a 
eutectic  particle  to  form.  The  transition  distance  of 
eutectic  particle  favorability  is  between  86  and  93  /tm, 
i.e.,  -1/3  of  the  PDAS  (248.4  fim).  This  plot  shows  that  it 


(a)  Original  image  (b)  Segmented  image  (c)  Euclidean  map 


Fig.  10 — (a)  A  1  mm  x  1  mm  subregion  from  Fig.  4  is  shown  along  with  two  corresponding  images  of  the  same  area:  ( h )  a  binary  image  with 
segmented  eutectic  particles  (white)  and  (c)  a  Euclidean  distance  map  from  the  dendrite  core  centroids,  where  lighter  intensity  refers  to  further 
distances  from  the  dendrite  cores. 
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Distance  from  Dendrite  Core  (pm) 

(a) 


(b) 


Fig.  11 — The  probability  of  a  eutectic  particle  as  a  function  of  the  distance  to  (a)  the  nearest  dendrite  core  or  ( b )  the  nearest  Voronoi  vertex. 
The  inset  schematic  shows  the  reference  point(s)  for  the  Euclidean  distance  in  each  plot. 


E» 


Fig.  12 — The  eutectic  particle  size  as  a  function  of  the  distance  to  (a)  the  nearest  dendrite  core  or  (b)  the  nearest  Voronoi  vertex.  The  distance 
for  each  particle  is  the  distance  for  the  particle  centroid.  The  50th-percentile  area,  A50  =  410  pm2,  refers  to  the  particle  size  where  50  pet  of  the 
eutectic  particle  area  lies  above/below  this  size. 


is  not  favorable  for  eutectic  particles  to  form  close  to  the 
primary  dendrite  core. 

Figure  11(b)  is  a  similar  plot  as  a  function  of  distance 
from  the  vertices  of  the  Voronoi  tessellation  (see 
schematic).  This  plot  was  generated  in  a  similar  manner 
to  Figure  11(a);  a  Euclidean  distance  map  was  first 
formed  from  the  Voronoi  vertices,  then  the  boundary 
pixels  within  100  /im  of  the  image  boundaries  were 
excluded,  etc.  There  is  an  increased  occurrence  of 
eutectic  particles  at  vertices,  regardless  of  their  distance 
from  the  dendrite  core.  This  observation  (along  with  the 
fact  that  the  probability  of  occurrence  is  higher  than  in 
Figure  11(a)  by  almost  2  pet)  suggests  that  solute  is 


forced  near  the  Voronoi  vertices,  thereby  increasing  the 
probability  of  eutectic  particle  occurrence.  The  transi¬ 
tion  distance  in  this  plot  is  between  67  and  90  ji m,  i.e.,  at 
-1/3  of  the  PDAS.  While  this  analysis  shows  the 
preferential  formation  of  eutectic  particles  based  on 
the  local  distances,  correlation  with  the  size  of  particles 
is  also  important. 

The  eutectic  particle  size  may  be  correlated  with  the 
distance  from  the  dendrite  cores  or  Voronoi  vertices  as 
well.  Figure  12  shows  the  eutectic  particle  size  as  a 
function  of  the  distance  from  the  nearest  dendrite  core 
and  the  nearest  Voronoi  vertice.  The  solid  line  shows  the 
50th-percentile  area,  ^450  =  410  /jm2,  which  refers  to  the 
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Fig.  13 — The  probability  of  a  eutectic  particle  of  a  certain  size  occurring  as  a  function  of  the  distance  to  (a)  the  nearest  dendrite  core  or  ( b )  the 
nearest  Voronoi  vertex.  Two  particle  sizes  are  considered:  particle  sizes  below  and  above  the  50th-percentile  area  A50.  The  solid  line  denotes  the 
distance  at  which  the  probability  functions  first  intercept,  indicating  a  transition  from  the  favorability  of  small  particles  to  large  particles  (a)  or 
vice  versa  (b). 


eutectic  particle  size  where  50  pet  of  the  eutectic  particle 
area  lies  above/below  this  size.  There  is  a  noticeable 
tendency  for  the  larger  particles  (A  >  A50)  to  form 
further  away  from  the  dendrite  core  and  closer  to  the 
Voronoi  vertices,  while  the  smaller  eutectic  particles 
(A  <  A50)  can  form  at  all  distances.  However,  it  is 
difficult  to  quantitatively  tell  from  the  following  plot 
what  the  preference  is  for  smaller  or  larger  particles  as  a 
function  of  distance.  Therefore,  to  further  quantify  this 
relationship  with  respect  to  the  size  of  the  particles,  the 
probability  associated  with  a  eutectic  particle  pixel 
belonging  to  either  a  small  or  large  particle  is  calculated 
in  Figure  13.  Interestingly,  in  Figure  13(a),  at  distances 
closer  to  the  dendrite  cores,  there  is  a  clear  preference 
for  smaller  particles  (A  <  A50 )  to  form  over  larger 
particles  (A  >  ^50).  At  a  distance  of  84.5 /im  (—1/3 
PDAS),  as  denoted  by  the  solid  line,  there  is  a 
crossover  in  the  probability  function  and  larger  parti¬ 
cles  are  statistically  favored  to  form  over  smaller 
particles.  In  the  case  of  distances  from  the  Voronoi 
vertices,  there  is  a  similar  behavior  except  that  larger 
particles  are  favored  at  smaller  distances  (closer  to 
Voronoi  vertices).  The  crossover  in  the  probability 
functions  occurs  at  79.3  /im  (-1/3  PDAS  again).  At 
distances  greater  than  this,  there  is  not  as  definitive  of  a 
trend  as  with  the  dendrite  cores,  i.e.,  in  some  cases, 
there  is  a  greater  probability  for  smaller  particles  to 
form  and,  in  some  cases,  for  larger  particles  to  form. 
Clearly,  the  distance  from  the  dendrite  cores  and, 
hence,  the  local  PDAS  affect  the  probabilities  of 
interdendritic  particles  to  form,  though.  In  a  similar 
manner,  it  is  anticipated  that  a  similar  relationship  may 
be  associated  with  shrinkage  porosity,  gas  porosity, 
and  other  interdendritic  defects. 


IV.  CONCLUSIONS 

Characterizing  the  PDAS  in  directionally  solidified 
microstructures  is  an  important  step  for  developing 
process-structure-property  relationships  by  enabling 
the  quantification  of  (i)  the  influence  of  processing  on 
microstructure  and  (ii)  the  influence  of  microstructure 
on  properties.  Thin-walled  directionally  solidified  struc¬ 
tures  ( e.g .,  a  turbine  blade)  require  new  approaches  for 
characterizing  the  dendrite  arm  spacing  and  the  micro¬ 
structure.  In  this  work,  we  utilized  a  new  Voronoi-based 
approach  for  spatial  point  pattern  analysis  that  was 
applied  to  an  experimental  dendritic  microstructure. 
This  technique  utilizes  a  Voronoi  tessellation  of  space 
surrounding  the  dendrite  cores  to  determine  nearest 
neighbors  and  the  local  PDAS.  In  addition,  we  com¬ 
pared  this  technique  to  a  recent  distance-based  tech¬ 
nique,  the  Warnken-Reed  method,  and  a  modification 
to  this  using  Voronoi  tessellations,  along  with  the 
minimal  spanning  tree  method.  Moreover,  a  convex 
hull-based  technique  was  used  to  include  edge  effects  for 
such  techniques,  which  can  be  important  for  thin 
specimens.  These  methods  were  used  to  quantify  the 
distribution  of  local  PDASs  as  well  as  their  spatial 
distribution  for  an  experimental  directionally  solidified 
superalloy  micrograph.  Last,  eutectic  particles  were 
segmented  to  correlate  distances  from  dendrite  cores 
and  Voronoi  vertices  to  the  occurrence  and  size  of  these 
interdendritic  features.  Interestingly,  with  respect  to  the 
distance  from  the  dendrite  core,  it  was  found  that  there 
is  a  greater  probability  of  occurrence  of  large  eutectic 
particles  (>410  /tm)  over  small  particles  at  dis¬ 
tances  greater  than  -1/3  of  the  bulk-measured  PDAS. 
This  systematic  study  of  the  different  techniques  for 
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quantifying  local  PDASs,  and  their  effect  on  micro¬ 
structure,  can  be  an  important  step  for  correlating 
with  both  processing  and  properties  in  single-crystal 
nickel-based  superalloys. 
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MARK  A.  TSCHOPP,  JONATHAN  D.  MILLER.  ANDREW  L.  OPPEDAL, 
and  KIRAN  N.  SOLANKI 

Microstructure  characterization  continues  to  play  an  important  bridge  to  understanding  why 
particular  processing  routes  or  parameters  affect  the  properties  of  materials.  This  statement 
certainly  holds  true  in  the  case  of  directionally  solidified  dendritic  microstructures,  where 
characterizing  the  primary  dendrite  arm  spacing  is  vital  to  developing  the  process-structure- 
property  relationships  that  can  lead  to  the  design  and  optimization  of  processing  routes  for 
defined  properties.  In  this  work,  four  series  of  simulations  were  used  to  examine  the  capability  of 
a  few  Voronoi-based  techniques  to  capture  local  microstructure  statistics  (primary  dendrite  arm 
spacing  and  coordination  number)  in  controlled  (synthetically  generated)  microstructures. 
These  simulations  used  both  cubic  and  hexagonal  microstructures  with  varying  degrees  of 
disorder  (noise)  to  study  the  effects  of  length  scale,  base  microstructure,  microstructure  vari¬ 
ability,  and  technique  parameters  on  the  local  PDAS  distribution,  local  coordination  number 
distribution,  bulk  PDAS,  and  bulk  coordination  number.  The  Voronoi  tesselation  technique 
with  a  polygon-side-length  criterion  correctly  characterized  the  known  synthetic  microstruc¬ 
tures.  By  systematically  studying  the  different  techniques  for  quantifying  local  primary  dendrite 
arm  spacings,  we  have  evaluated  their  capability  to  capture  this  important  microstructure 
feature  in  different  dendritic  microstructures,  which  can  be  an  important  step  for  experimentally 
correlating  with  both  processing  and  properties  in  single  crystal  nickel-based  superalloys. 
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I.  INTRODUCTION 

MlCROSTRUCTURE  imaging  and  characterization 
are  an  important  bridge  to  understanding  why  par¬ 
ticular  processing  routes  or  parameters  affect  the 
properties  of  materials.  Developing  an  enhanced  under¬ 
standing  of  mechanical  behavior  in  materials  relies  upon 
sufficiently  characterizing  microstructure  details  at  the 
relevant  length  scales  that  contribute  to  this  behavior. 
Moreover,  to  truly  enhance  the  predictive  capability  of 
processing-structure-property  models  that  aim  to 
improve  material  performance  requires  a  quantitative 
stereological  description  of  the  relevant  microstructure 
features  and,  thereby,  the  material  itself.  Predictive 
models  that  effectively  capture  the  linkage  between 
processing  and  properties  (through  microstructure)  can 
be  utilized  within  an  integrated  computational  materials 
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engineering  approach  to  design  materials  and  accelerate 
their  insertion  into  application. 

As  crystallographic  texture  and  grain  boundary  prop¬ 
erties  are  instrumental  to  the  properties  and  perfor¬ 
mance  of  some  materials,  processing  routes  that  produce 
directionally  solidified  dendritic  microstructures  have 
emerged  as  a  way  to  remove  adverse  effects  of  grain 
boundaries  that  are  observed  in  some  materials  appli¬ 
cations.  One  such  application  is  in  turbine  blades  within 
the  high  temperature  section  of  the  modern  turbine 
engine,  where  single  crystal  nickel-based  superalloys  are 
widely  used.11’21  In  single  crystal  nickel-based  superal¬ 
loys,  there  are  a  number  of  length  scales  of  microstruc- 
ture  that  contribute  to  mechanical  behavior,  ranging 
from  the  /  precipitates  to  pores  and  eutectic  particles  to 
the  dendrites  themselves.  At  the  largest  microstructure 
length  scale  in  directionally  solidified  single  crystal 
microstructures,  the  features  of  interest  are  the  den¬ 
drites;  many  features  at  lower  length  scales  (e.g.,  eutectic 
particles,  precipitates)  or  at  similar  scales  (e.g.,  porosity, 
freckle  defects)  are  strongly  associated  with  the  dendrite 
arm  spacing  and  morphology. 13-71  The  solidification 
morphology  associated  with  dendrite  arm  spacing  has 
been  described  in  the  literature.18'91 

Historically,  the  primary  dendrite  arm  spacing 
(PDAS)  has  been  found  to  correlate  with  processing 
(e.g.,  solidification/cooling  rate)17' 10-1 71  as  well  as  with 
properties  (e.g.,  creep  strength,  fatigue  proper¬ 
ties).16-16'1  S1  Lamm  and  Singer161  produced  single  crystal 
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nickel-based  microstructures  (PWA  1483)  with  a  varied 
range  of  different  dendrite  arm  spacings  (250  to  600 
/im).  In  their  study,  they  found  that  decreasing  the 
mean  dendrite  arm  spacing  was  associated  with  an 
increased  high-cycle  fatigue  life,  because  the  fatigue 
cracks  tended  to  originate  at  shrinkage  porosity;  the 
largest  shrinkage  pores  correlated  with  a  large  PDAS. 
Similar  observations  of  a  correlation  between  dendrite 
arm  spacings  and  microporosity  or  solute  segregation  has 
been  observed  by  researchers  as  well.17,151''  On  the 
processing  side,  manufacturing  turbine  blades  via  the 
liquid  metal  cooling  (LMC)  process14,14’19,201  enables  a 
greater  degree  of  control  over  the  size  of  the  primary 
dendrite  arm  spacing  because  of  the  higher  thermal 
gradient  associated  with  this  process.  For  example, 
Franke  et  alp 11  summarizes  several  recent  studies  com¬ 
paring  the  LMC  process  and  the  conventional  high  rate 
solidification  process.  In  each  case  considered,  the  LMC 
process  resulted  in  a  smaller  primary  dendrite  arm 
spacing  and  an  improved  withdrawal  speed.  Additionally, 
some  researchers  are  exploring  the  use  of  magnetic  fields 
to  complement  directional  solidification  processes  and 
further  reduce  the  primary  dendritic  arm  spacing.  122-241 
The  use  of  various  processing  techniques  such  as  LMC 
may  allow  thermal  gradients,  and  hence  primary  dendrite 
arm  spacing,  to  be  tailored  towards  a  (more  desirable) 
lower  mean  primary  dendrite  arm  spacing  and  more 
homogeneous  primary  dendrite  arm  spacing.1141 

The  identification  of  the  primary  dendrites  in  these 
studies  is  required  to  quantify  the  effect  of  processing/ 
alloying  on  the  microstructure.  In  many  experimental 
studies  of  dendrite  arm  spacing,  the  dendrite  core  is 
typically  identified  manually  in  the  experimental  micro¬ 
graphs.  In  light  of  some  of  the  larger  dendrite  datasets 
and  images,  there  are  some  automated  image  processing 
methods  that  have  been  investigated  to  identify  dendrite 
core  locations  by  exploiting  their  fourfold  symmetry.1251 
Some  researchers1261  have  highlighted  some  of  the 
difficulties  in  identifying  dendrite  arm  spacing  from 
microstructure  and  have  suggested  modified  alternatives 
based  on  dendrite  trunk  diameters,1271  which  may  also  be 
useful  for  correlating  with  solidification  processing 
parameters.  Hence,  for  understanding  the  role  of  den¬ 
drite  arm  spacing  in  future  alloys  and  processing  routes, 
a  confluence  of  automated  experimental  techniques  for 
attaining  microstructures,  automated  techniques  for 
identifying  dendrites,  and  automated  (and  informed) 
characterization  techniques  is  required. 

The  traditional  approach  for  measuring  primary 
dendrite  arm  spacing  in  single  crystal  metals,  whereby 
the  number  of  dendrite  cores  in  a  specified  area  is  related 
to  the  dendrite  arm  spacing110,28’291  is  given  by: 

X  =  cJ±  [1] 

V  n 

where  X  is  primary  dendrite  arm  spacing,  A  is  the  area 
analyzed,  n  is  the  number  of  dendrites,  and  c  is  a 
coefficient  that  depends  on  the  microstructure.  McCart¬ 
ney  and  Hunt1101  showed  that  c  =  0.5  for  a  random 
array  of  points,  c  =  1  for  a  square  array  of  points,  and 
c  =  1.075  for  a  hexagonal  array  of  points;  they  had  to 
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apply  a  correction  for  the  bulk  dendrite  arm  spacing  X  as 
processing  conditions  caused  a  change  in  the  local 
environment  of  the  dendrites,  i.e.,  a  transition  from 
cellular  growth  to  dendrite  growth.110,1 1,301  However, 
this  approach  is  insufficient  for  capturing  local  arm 
spacings  or  the  dendrite  arm  spacing  distribution,  and 
may  provide  problems  with  complex  geometries  such  as 
turbine  blades.  In  fact,  part  of  the  motivation  for 
quantifying  the  local  PDAS  is  that  a  narrow  distribution 
(i.e.,  low  standard  deviation)  of  local  PDAS  values  may 
result  in  a  more  homogeneous  distribution  of  interden- 
dritic  microstructure  features  and,  more  importantly,  a 
narrow  distribution  of  mechanical  properties.1311  Struc¬ 
tural  inhomogeneities  in  single  crystal  components, 
which  may  cause  wider  property  variations,  have  been 
correlated  to  asymmetric  heat  flux  and  transients  in 
solidifications.1321 

Previous  work  examined  some  recent  approaches,  as 
well  as  some  modified  versions  of  these  approaches,  for 
characterizing  the  local  dendrite  arm  spacing  within 
experimental  dendritic  microstructures.1  1  This  work 
used  an  experimental  dendritic  microstructure  along 
with  three  different  techniques  that  are  based  on  the 
nearest  neighbor  spacing  and/or  a  Voronoi  tesselation 
of  the  dendrite  cores.  Figure  1  shows  (a)  an  example  of  a 
directionally  solidified  dendritic  microstructure1341  along 
with  contour  plots  showing  the  evolution  of  the  local 
PDAS  as  a  function  of  dendrite  location  for  four 
different  techniques  (to  be  discussed  in  Section  II  B). 
This  study  found  that  augmenting  existing  techniques 
with  Voronoi  tesselations  resulted  in  a  more  physical 
description  of  the  local  dendrite  arm  spacing  and 
resulted  in  a  mean  PDAS  that  compared  favorably  to 
the  traditional  bulk  characterization  technique  (Eq.  [1]). 
However,  several  questions  remained.  For  example, 
while  comparing  the  mean  of  the  local  PDAS  with  the 
bulk  PDAS  calculated  from  Eq.  [1]  is  a  valid  method  for 
examining  the  validity  of  the  local  characterization 
techniques,  an  experimental  microstructure  contains 
some  degree  of  disorder.  Hence,  it  is  not  clear  whether 
these  techniques  are  accurately  capturing  both  the  local 
PDAS  and  the  local  coordination  number,  since  an 
experimental  microstructure  may  contain  regions  of 
fourfold  coordination  and  sixfold  coordination.  Since 
the  coordination  number  may  bias  the  local  PDAS 
calculation,  it  is  important  to  understand  the  extent  to 
which  these  techniques  capture  both  local  measures. 

Therefore,  the  research  objective  of  the  present  work 
is  to  understand  the  ability  of  these  Voronoi-based 
techniques  to  capture  local  microstructure  statistics 
(primary  dendrite  arm  spacing  and  coordination  num¬ 
ber)  in  controlled  (synthetically  generated)  microstruc¬ 
tures.  To  meet  this  objective,  cubic  and  hexagonal 
synthetic  microstructures  with  varying  degrees  of  disor¬ 
der  (noise)  were  generated  to  study  the  effects  of  length 
scale,  base  microstructure,  microstructure  variability, 
and  technique  parameters  on  the  local  PDAS  distribu¬ 
tion,  local  coordination  number  distribution,  bulk 
PDAS,  and  bulk  coordination  number  (calculated  as 
the  mean  of  the  local  statistics). 

The  paper  is  organized  as  follows.  Section  II  defines 
the  methodology  used  herein,  including  the  methods 
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Fig.  1 — (a)  Example  of  dendritic  structure  normal  to  the  withdrawal  direction  in  a  directionally  solidified  single  crystal  nickel-based  superalloy 
cast  using  the  liquid  metal  cooling  technique)3^  The  dots  denote  the  dendrite  cores,  where  the  white  dots  indicate  dendrite  cores  at  the  edge  of 
the  field  of  view.  The  local  primary  dendrite  arm  spacing  (/an)  was  calculated  using  four  different  techniques:^  (h)  Voronoi  tesselation  with 
t/crit  =  0.12,  (c)  Voronoi  tesselation  with  t/crjt  =  0.00,  (d)  Warnken-Reed  technique  with  a  =  2.0,  and  (e)  Voronoi  Warnken-Reed  with  a  —  2.0. 


used  to  generate  the  synthetic  microstructures,  the 
different  techniques  used  to  quantify  the  local  PDAS, 
and  the  simulation  test  matrix  that  was  performed. 
Section  III  presents  the  results,  analysis,  and  discussion 
based  on  the  various  sets  of  simulations  that  were 
performed.  This  includes  studies  on  the  effect  of  the 
microstructure  length  scale;  the  effect  of  microstructure 
variability  and  technique  parameters  on  the  bulk 
responses;  the  effect  of  microstructure  variability  and 
technique  parameters  on  the  local  distributions;  and  the 
correlation  between  different  local  measures.  Section  IV 
presents  the  conclusions.  The  significance  of  the  present 
work  is  that  it  aims  to  answer  several  questions 
pertaining  to  how  to  measure  the  local  primary  dendritic 
arm  spacing  in  experimental  dendritic  microstructures. 


H.  METHODOLOGY 

A.  Synthetic  Microstructure  Generation 

Synthetic  microstructures  were  used  to  assess  the 
ability  of  the  different  methods  examined  within  to 
capture  perturbations  from  two  symmetric  arrange¬ 
ments:  the  4-nearest  neighbor  cubic  cell  and  the 
6-nearest  neighbor  hexagonal  cell.  In  addition  to  the 
perturbations  from  the  ideal  cubic  and  hexagonal 
configurations,  a  microstructure  with  complete  random¬ 
ness  of  points  was  generated  using  random  number 
generators  to  create  N  independent  points  with  coordi¬ 


nates  Rx  and  Rv  for  their  respective  x  and  y  coordinates 
within  the  appropriate  bounds.  Figures  2(a)  through  (c) 
illustrates  the  three  different  types  of  patterns  used  in  the 
present  study.  The  random  pattern  was  added  to 
quantify  when  the  perturbations  for  a  specific  pattern 
became  indistinguishable  front  a  random  set  of  points. 

The  degree  of  deviation  from  the  ideal  symmetric 
configurations  was  then  systematically  introduced.  Pre¬ 
vious  methods  have  used  a  uniform  or  Gaussian 
distribution  of  noise  in  both  the  x  and  y  directions  for 
all  points  or  have  utilized  a  combination  of  the  fraction 
of  points  perturbed  as  well  as  the  distance  shifted.  Here, 
we  have  utilized  a  uniform  distribution  of  noise  that  is  a 
fraction  of  the  nearest  neighbor  distance  oq.  Figures  2(d) 
through  (f)  shows  three  different  noise  levels  for  a  10  x 
10  cubic  pattern,  whereby  the  lines  illustrate  the  degree 
of  deviation  from  the  starting  base  lattice.  As  noise  is 
increased,  the  starting  lattice  becomes  less  recognizable 
as  it  moves  towards  a  purely  random  configuration  of 
points.  However,  the  ability  to  characterize  the  local 
environment  of  a  point  in  the  presence  of  these 
perturbations  and  relate  this  to  a  base  microstructure 
lattice  presents  a  number  of  challenges.  The  subsequent 
subsection  will  discuss  the  techniques  that  are  used 
herein  to  characterize  such  microstructures. 

Dendrites  near  the  free  surface  complicate  the  analysis 
of  the  PDAS  statistics.  These  have  been  addressed 
previously  using  a  convex  hull  and  Boolean  operations 
to  enable  computing  of  localized  statistics  using  Voronoi 
tesselations.  Figure  2(g)  shows  both  interior  points 
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Fig.  2 — Three  different  types  of  patterns  were  used  in  the  present  study:  (a)  cubic  pattern,  ( b )  hexagonal  pattern,  and  (c)  random  pattern.  Three 
different  microstructures  show  the  effect  of  the  noise  level  of  ( d)  0.20ao,  (^)  0.40ao,  and  (/)  0.60«o-  As  an  illustration  of  the  deviation  from  the 
perfect  cubic  array  of  points,  the  original  coordinates  of  the  cubic  pattern  (i.e.,  Fig.  2(a))  is  shown  as  a  black  dot  with  a  line  connecting  it  to  the 
deviated  coordinate  (large  blue  dots).  For  this  particular  10  x  10  pattern,  as  the  noise  level  is  increased,  the  same  points  are  shifted  in  the  same 
direction;  the  only  change  is  in  the  magnitude  of  the  shift,  (g)  To  eliminate  edge  effects,  statistics  from  points  on  the  exterior  of  the  area  are  ex¬ 
cluded.  However,  these  points  are  still  used  to  evaluate  statistics  in  points  on  the  interior  region. 


(blue)  and  edge  points  (red)  for  the  large  microstruc¬ 
tures,  whereby  the  edge  points  are  selected  using  a  small 
distance  (1.5«o)  front  edge  of  the  microstructure  pattern. 
The  present  analysis  includes  these  edge  points  only  for 
calculating  the  localized  statistics  of  interior  points,  but 
does  not  use  their  localized  statistics  for  subsequent 
mean  statistics,  heretofore  also  referred  to  as  the  bulk 
statistics.  In  this  sense,  the  present  paper  focuses  on 
locally  measured  statistics  and  their  mean  responses 
(i.e.,  the  bulk  PDAS  and  bulk  coordination  number) 
within  a  bulk  microstructure. 

Last,  it  should  be  mentioned  that  there  are  no  rules 
during  the  generation  step  that  prohibit  these 
microstructures  front  having  two  neighboring  points  that 
are  in  close  proximity  to  one  another  (e.g.,  the  random 
microstructure  in  Figure  2(c)).  In  an  experimental  mi¬ 
crostructure,  the  primary  dendrite  cores  are  separated  by 
a  critical  minimum  distance  that  should  scale  with  both 
the  size  of  the  dendrite  stalk  and  the  distance  associated 
with  the  transition  of  tertiary  to  primary  dendrite  arms. 
For  example,  our  previous  study^  shows  that  the 
dendrite  stalk  diameter  is  approximately  20  pet  of  the 
bulk  PDAS  (i.e.,  stalk  diameter  is  twice  the  maximum 
distance  at  which  no  eutectic  particles  were  observed). 
This  prior  result  suggests  that  noise  levels  of  0.40uo  are 
the  highest  noise  level  that  might  be  observed  in  a  real 
microstructure,  because  any  higher  noise  level  could 
produce  neighboring  points  that  lie  within  0.20ao  distance 
of  each  other.  However,  this  estimation  of  the  trunk 
diameter  is  only  an  estimation — other  studies  have  shown 
that  the  trunk  diameter  itself  may  be  a  useful  criterion 
for  understanding  directional  solidification  and  its  impact 
on  microstructureJ271  While  not  examined  herein, 
microstructures  built  using  a  shell  model  or  random 
sequential  adsorption  algorithms  (e.g.,  Reference  35) 
could  allow  researchers  to  understand  how  this  critical 


minimum  distance  affects  the  local  PDAS  results.  In 
general,  though,  the  higher  degrees  of  noise  in  the 
synthetic  microstructures  allow  us  (1)  to  capture  and 
understand  the  trends  in  mean  PDAS  behavior,  and 
(2)  test  the  capability  of  these  different  algorithms  to 
handle  some  potentially  extreme  scenarios  (of  2  or  more 
points  lying  close  to  one  another).  Furthermore,  the 
results  show  that  the  random  microstructures  primarily 
comprised  a  mixture  of  4-,  5-,  and  6-coordinated  points, 
in  agreement  with  experimental  microstructures. 

B.  Local  Characterization  Techniques 

The  approach  utilized  herein  to  measure  the  local 
dendrite  arm  spacing  is  based  on  a  Voronoi  tesselation 
of  the  spatial  array  of  dendrite  cores.  The  following 
analysis  techniques  were  implemented  in  MATLAB 
R2014  (The  MathWorks,  Inc.).  To  illustrate  how  the 
present  method  works  and  differs  from  some  other 
published  methods,  a  synthetic  5x5  cubic  pattern  of 
points  is  generated  with  varying  degrees  of  noise 
(100  pet  noise  fraction  and  noise  fraction  as  in  Refer¬ 
ences  36,  37),  as  shown  in  Figure  3.  In  Figure  3,  the 
different  techniques  for  the  same  microstructure  are 
shown  from  left  to  right  (labels  at  the  bottom),  while  the 
noise  level  is  increased  from  top  to  bottom  to  illustrate 
the  differences.  The  techniques  presented  in  Figure  3 
were  first  used  together  in  Tschopp  et  alP 3'  to  charac¬ 
terize  the  local  spacing  for  an  experimental  microstruc¬ 
ture  of  a  directionally  solidified  single  crystal  nickel- 
based  superalloy.  In  contrast  to  that  work,  the  following 
work  will  examine  the  sensitivity  and  applicability  of  a 
few  of  these  techniques  based  on  thousands  of  synthetic 
microstructure  instantiations. 

Figure  3  shows  several  methods  for  characterizing  the 
local  microstructure  environment.  One  such  method  for 
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Fig.  3 — The  difference  between  various  methods  for  defining  the  nearest  neighbors  (red  dots)  and  their  spacing  for  a  single  point  (large  black 
dot)  with  a  distorted  local  environment.  Images  in  (a)  through  (/)  are  (a)  the  initial  pattern,  (b)  the  Warnken-Reed  method  with  a  =  1.5  and 
^initial  =  3,  (c)  the  Voronoi  tesselation  diagram  for  the  points,  (d)  the  Voronoi-modified  Warnken-Reed  method  with  a=  1.5  and  initial  =  3, 
(e)  the  Voronoi  method  (f/crjt  =  0.0),  and  (/)  the  modified  tesselation-based  technique  with  dcrn  =  0.10. 


measuring  the  local  dendrite  arm  spacing  is  the 
Warnken-Reed  methodi36-37!  The  Warnken-Reed 
method  calculates  the  dendrite  arm  spacing  for  a  single 
point  (black  dot)  by  starting  with  an  initial  number  of 
nearest  neighbors,  initial  (often  =  3,  i.e.,  three 

closest  neighbors),  and  then  iteratively  adding  potential 
nearest  neighbors  that  are  within  a  cutoff  distance 
defined  by  the  already-added  nearest  neighbors.  The 
inner  circle  in  the  second  and  fourth  column  of  Figure  3 


represents  the  average  spacing,  davg,  of  these  neighbors 
and  the  outer  circle  represents  the  cutoff  for  adding  the 
next  neighbor,  rfavg  +  arfstd,  where  dst a  is  the  standard 
deviation  of  the  nearest  neighbor  spacings  and  a  is  a 
parameter  that  is  typically  between  1  and  2.  Neighbors 
continue  to  be  added  until  the  cutoff  does  not  include 
any  new  neighbors.  The  local  coordination  number  and 
dendrite  arm  spacing  is  calculated  from  the  neighbors 
added  (shown  as  red  dots  in  Figure  3).  However,  if  the 


4614— VOLUME  46A,  OCTOBER  2015 


METALLURGICAL  AND  MATERIALS  TRANSACTIONS  A 


41 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


standard  deviation  of  the  distances  of  the  nearest 
neighbors  dsU j  or  the  parameter  a  is  large,  this  technique 
can  continue  to  add  nearest  neighbors  beyond  the  first 
nearest  neighbors.  Notice  that  the  distance  between  the 
inner  and  outer  circle  tends  to  increase  with  increasing 
noise  due  to  the  increased  r/std  of  the  nearest  neighbors. 
As  discussed  in  Tschopp  el  a/.,'33'  clearly  a  method  for 
restricting  the  number  of  nearest  neighbors  using  such  a 
technique  is  necessary. 

Another  method  of  identifying  the  potential  of  first 
nearest  neighbors  is  to  perform  a  Voronoi  tesselation  of 
the  space  surrounding  the  points,  as  shown  in  the  third 
column  of  Figure  3.  The  polygon  edges  are  equidistant 
between  the  points  contained  in  the  two  adjacent 
polygons  and  the  triple  points  (merging  of  three  lines) 
are  equidistant  between  the  points  contained  in  the  three 
adjacent  polygons.  Therefore,  the  first  nearest  neighbors 
(FNNs,  shown  as  open  circles  in  Figure  3)  correspond  to 
the  edges  of  the  central  polygon  (that  contains  the  black 
dot).  This  subset  of  points  is  the  maximum  number  of 
nearest  neighbors  that  the  central  point  can  have. 

In  this  manner,  several  techniques  have  been  identi¬ 
fied  to  quantify  a  local  dendrite  arm  spacing  based  on 
the  Voronoi-identified  FNNs.1331  For  example,  the 
Voronoi  Warnken-Reed  (VWR)  method  (Figure  3) 
only  includes  the  Voronoi  FNNs  as  potential  nearest 
neighbors  and  cannot  expand  beyond  these,  alleviating  a 
potential  problem  of  selecting  second-nearest  neighbors 
or  greater.  Another  method  using  the  Voronoi  FNNs  is 
to  consider  all  of  these  potential  nearest  neighbors  as 
nearest  neighbors  (Figure  3),  as  in  Brundidge  et  a/.'15'. 
Unfortunately,  this  approach  is  sensitive  to  small 
perturbations  in  the  spatial  positions  of  the  neighbors. 
Notice  that  a  small  degree  of  noise  added  to  the  cubic 
lattice  (in  row  2)  results  in  several  more  nearest 
neighbors  being  identified.  As  the  corresponding  edge 
of  the  polygon  for  these  45  deg  points  is  small,  a  slight 
movement  away  from  the  central  point  would  cause  this 
point  to  no  longer  share  an  edge  with  the  polygon 
containing  the  black  dot;  in  this  scenario,  the  two 
adjacent  polygons  on  either  side  would  effectively 
“pinch"  off  this  neighbor.  This  scenario,  however,  has 
a  physical  basis  in  solidification  microstructures  as  these 
can  be  envisioned  as  two  dendrite  cores  that  compete 
with  the  central  core,  and  the  45  deg  cores  have  a  much 
less  prominent  effect  on  the  central  core. 

The  last  method  utilizes  a  criterion  based  on  the  edge 
lengths  of  the  Voronoi  polygon.133'  In  Figure  3,  those 
neighbors  with  edge  lengths  less  than  a  critical  fraction. 
dCI it,  of  the  total  polygon  perimeter  are  excluded  as  nearest 
neighbors  (e.g.,  10  pet  in  Figure  3).  In  the  present  study, 
the  local  dendrite  arm  spacing  statistics  are  evaluated 
using  two  of  these  four  techniques:  VWR  and  the  Voronoi 
technique  with  (r/crjt  >  0)  a  line  length  threshold. 

With  respect  to  the  addition  of  noise  to  the  5x5 
cubic  pattern,  these  methods  handle  the  increased  levels 
of  noise  very  differently.  First,  the  distance-based 
Warnken-Reed  method  does  not  incorporate  enough 
neighbors  at  low  noise  levels  and  incorporates  many 
more  neighbors  at  higher  noise  levels  than  the  Voronoi- 
based  techniques  (rightmost  two  columns).  The  use  of 
the  Voronoi-identified  nearest  neighbors  to  limit  the 


number  of  neighbors  selected  does  seem  to  improve  the 
technique  in  a  physical  manner  for  higher  noise  levels. 
For  several  noise  levels,  the  Voronoi  method  with 
r/crit  =  0.0  does  clearly  overestimate  the  number  of 
nearest  neighbors,  while  the  four  nearest  neighbors 
identified  through  dCI jt  =  0.1  (Figure  3)  perhaps  offers  a 
better  approximation  of  the  number  of  nearest  neigh¬ 
bors.  Interestingly,  the  Voronoi  method  with  dCT jt  =  0.1 
does  identify  four  nearest  neighbors  for  all  noise  levels. 
Previous  work  has  also  shown  that  while  these  methods 
can  identify  the  same  number  of  neighbors,  they  can  in 
fact  be  different  neighbors.'33'  This  can  be  due  to  the 
Warnken-Reed  method  being  a  distance-based  method, 
and  identifying  the  four  closest  neighbors,  while  the 
modified  Voronoi  technique  is  based  on  the  edge  lengths 
of  the  Voronoi  polygon,  and  hence  utilizes  this  to 
identify  nearest  neighbors  (which  may  not  be  the  closest 
neighbors). 

The  motivation  for  examining  the  sensitivity  and 
capability  of  these  local  PDAS  metrics  is  that  the 
traditional  PDAS  metric  does  not  consider  the  order  or 
disorder  of  the  dendrites  within  the  microstructure. 
Figure  3  illustrates  why  a  local  metric  for  PDAS  may  be 
needed.  For  the  field  of  view  given  in  Figures  2  and  3, 
the  bulk  PDAS  metric  would  be  the  same  since  the 
number  of  dendrites  n  and  the  area  A  are  equal  (see 
Eq.  [1]).  However,  the  increasing  disorder  of  the  den¬ 
dritic  structure  in  the  case  of  Figure  3  may  yield  (i)  a 
more  uneven  distribution  of  solute  elements,  (ii)  the 
formation  of  second  phase  particles,  (iii)  the  formation 
of  gas  or  shrinkage  porosity,  or  (iv)  the  lateral  growth  of 
secondary  dendrite  arms.  Hence,  in  addition  to  the  bulk 
PDAS  values,  understanding  how  processing  conditions 
may  impact  the  disorder  of  the  dendritic  structure  may 
be  important  for  understanding  the  properties  of  direc¬ 
tionally  solidified  alloys. 


C.  Simulation  Test  Matrix 

The  ability  to  rapidly  generate  synthetic  microstruc¬ 
tures  of  known  lattice  spacing,  noise  fraction,  and 
coordination  (from  the  type  of  pattern  utilized)  allows 
for  a  controlled  understanding  of  how  well  these 
computational  techniques  can  estimate  the  local  coor¬ 
dination  numbers  and  local  PDAS.  First,  several  studies 
were  explored  to  understand  how  these  different  factors 
affect  the  bulk  dendrite  responses  (PDAS  and  coordi¬ 
nation  number)  based  on  the  mean  value  of  the  locally 
calculated  statistics;  in  this  manner,  these  mean  respons¬ 
es  can  be  directly  compared  with  the  known  bulk  PDAS 
and  the  coordination  number  (for  a  zero-noise  mi¬ 
crostructure).  Next,  a  series  of  simulations  were  per¬ 
formed  to  examine  the  local  response  and  how  this 
changes  as  a  function  of  these  different  factors.  In  doing 
so,  the  present  work  aims  to  answer  several  questions 
pertaining  to  how  to  measure  the  local  primary  dendritic 
arm  spacing  in  experimental  microstructures,  as  posed 
in  the  following  4  studies. 

1.  What  is  the  effect  of  the  microstructure  length  scale 
on  bulk  (or  mean)  PDAS  measurements?  In  par¬ 
ticular,  the  goal  of  this  analysis  is  to  understand 
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how  changing  the  number  of  dendrites  (points)  cap¬ 
tured  within  an  image  (field  of  view)  relates  to  the 
overall  dendrite  arm  spacing.  This  set  of  simula¬ 
tions  varied  several  parameters  to  understand  this 
variability:  the  length  scale  of  the  field  of  view  (i.e., 
the  number  of  points/dendrites),  the  type  of  pattern 
used  (Figure  2),  and  the  type  of  method  used 
(VWR,  Voronoi).  Twenty  microstructure  instan¬ 
tiations  were  used  at  a  noise  fraction  of  0.40ao 
(a0  =  1). 

2.  What  is  the  effect  of  the  parameters  for  the  different 
methods  (VWR,  Voronoi)  on  the  bulk  PDAS  mea¬ 
surements?  How  does  the  variability  in  the  microstruc¬ 
ture  affect  these  measurement  techniques?  The  goal  of 
this  analysis  is  to  understand  how  a  and  c/crjt  for  the 
VWR  and  Voronoi  methods  affects  the  measured 
PDAS  and  coordination  number  for  different 
microstructures.  Hence,  a  was  varied  from  1.0  to  2.0 
in  increments  of  0.1;  similarly,  t/crjt  was  varied  from 
0.0  to  0.12  in  increments  of  0.02.  The  microstructure 
may  also  influence  the  sensitivity  of  these  parameters, 
so  both  cubic  and  hexagonal  microstructures  were 
generated  with  6  different  noise  levels:  0.05ao,  O.lOao, 
0.20ao,  0.30tfo,  0.40cro.  and  0.50ao- 

3.  What  is  the  effect  of  the  parameters  (a  and  r/crjt)  for 
the  different  methods  (VWR  and  Voronoi  Tessela- 
tion  techniques,  respectively)  on  the  local  PDAS 
and  the  local  coordination  number  distributions? 
How  is  this  affected  by  the  base  configuration  and 
noise  of  the  microstructure?  The  goal  here  is  to 
delve  further  than  the  mean  statistics  and  start  to 
understand  how  these  different  techniques  and  their 
parameters  affect  the  distributions  of  the  local 
statistics.  Hence,  the  same  simulations  as  described 
in  study  #2  will  be  used,  but  further  analysis  will  be 
performed  pertaining  to  the  local  distributions  and 
how  they  evolve. 

4.  What  is  the  correlation  between  local  PDAS  calcu¬ 
lated  at  two  different  parameters  for  the  same 
technique  and  same  microstructures?  What  is  the 
correlation  between  the  local  PDAS  and  the  local 
coordination  number  for  microstructures  with  dif¬ 
ferent  noise  levels?  How  is  the  bulk  PDAS  linked  to 
the  coordination  number  and  noise  level?  The  goal 
of  this  analysis  is  to  understand  the  changes  in  the¬ 
se  local  statistics  with  the  parameters  for  these  tech¬ 
niques.  In  other  words,  how  sensitive  is  any  local 
PDAS  to  the  parameter  chosen?  Also,  does  the 
coordination  number  directly  (or  linearly)  correlate 
with  the  calculated  PDAS?  This  study  will  investi¬ 
gate  these  relationships. 

In  short,  the  following  set  of  simulations  are  per¬ 
formed  on  different  synthetically  generated  microstruc¬ 
tures  with  varying  degrees  of  disorder  to  better 
understand  the  capability  of  these  sorts  of  local 
microstructure  characterization  techniques  as  it  applies 
to  the  problem  of  measuring  the  local  primary  dendrite 
arm  spacing  in  single  crystal  microstructures.  That  being 
said,  these  techniques  are  not  restricted  to  this  applica¬ 
tion,  but  can  also  be  applied  to  characterizing  local 
statistics  in  other  systems  and  problems  as  well. 
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III.  RESULTS  AND  DISCUSSION 

A.  Microstructural  Length  Scale  Effects  on  PDAS 

The  PDAS  was  quantified  for  multiple  microstructure 
instantiations  to  understand  the  effect  of  length  scales 
on  two  measurement  techniques.  Twenty  different  mi¬ 
crostructures  were  generated  for  six  different  mi¬ 
crostructure  lengths:  10  (10  x  10,  i.e.,  100  points),  20, 
30,  40,  50,  and  80.  In  this  example,  the  spacing  between 
nearest  neighbor  points/dendrites  was  chosen  to  be  1, 
such  that  the  PDAS  statistics  would  have  a  mean  value 
of  1  for  a  noise  fraction  of  0.  Recall  that  to  eliminate 
edge  effects,  statistics  from  points  on  the  exterior  of  the 
microstructure  are  excluded  (within  1.5  units  of  the 
edges).  Three  patterns  (cubic,  hexagonal,  and  random) 
were  used  to  assess  how  the  type  of  base  pattern  affects 
PDAS  convergence.  The  noise  fraction  was  chosen  to  be 
0.40ao,  i.e.,  large  enough  to  add  a  significant  amount  of 
variability,  but  not  so  large  as  to  potentially  cause  points 
to  overlap  (i.e.,  0.50ao  and  above).  Two  different 
techniques — VWR  and  Voronoi — were  used  to  assess 
the  convergence  of  the  PDAS  statistics  with  length  scale. 
The  VWR  method  used  a  =  2  and  the  Voronoi  tech¬ 
nique  used  dc Hi  =  0.12.  While  these  techniques  can  be 
used  to  calculate  the  local  PDAS  and  coordination 
number,  in  this  case  these  techniques  are  used  to 
calculate  a  bulk  PDAS  value,  which  is  the  mean  of  the 
local  PDAS  values. 

The  evolution  of  the  bulk  PDAS  and  the  standard 
deviation  of  the  bulk  PDAS  as  a  function  of  the  length 
scale  is  plotted  in  Figure  4.  In  Figure  4(a),  there  are 
several  features  to  notice.  First,  a  dotted  line  is  drawn  at 
a  PDAS  of  1  to  indicate  the  expected  PDAS  value  with  a 
noise  fraction  of  0.  Second,  each  data  point  is  the  mean 
value  of  the  20  bulk  PDAS  quantified  for  each  mi¬ 
crostructure  instantiation.  Similarly,  the  errors  bars  are 
equivalent  to  one  standard  deviation  of  those  20 
microstructure  instantiations.  The  error  bars  are  offset 
slightly  for  each  of  the  different  prescribed  length  scales 
to  avoid  the  significantly  overlap  and  aid  in  visualiza¬ 
tion.  Third,  the  different  symbols  indicate  the  different 
pattern  types  used  and  the  different  colors  delineate  the 
different  techniques  utilized. 

There  are  several  trends  evident  from  Figure  4(a). 
First,  as  the  length  scale  (and  size)  of  the  microstructure 
increases,  the  bulk  PDAS  values  converge  quickly,  such 
that  the  bulk  PDAS  value  is  very  similar  in  the  20  x  20 
microstructures  to  that  in  the  80  x  80  microstructures, 
as  would  be  expected.  However,  the  error  bars  denote 
that  there  is  significant  variation  in  the  individual  bulk 
PDAS  values,  which  decreases  as  the  length  scale  of  the 
microstructure  field  of  view  is  increased.  In  fact. 
Figure  4(b)  plots  the  standard  deviation  of  the  bulk 
PDAS  values  as  a  function  of  the  microstructure  length 
scale  on  a  log-log  scale  and  shows  a  similar  decrease  for 
all  synthetic  microstructures;  a  similar  decreasing  trend 
was  found  for  the  standard  deviation  of  the  bulk 
coordination  number  as  well.  In  addition  to  the  trends 
with  respect  to  the  length  scales,  which  is  somewhat 
independent  of  the  type  of  microstructure  and  measure¬ 
ment  technique,  there  are  trends  with  respect  to  the 


METALLURGICAL  AND  MATERIALS  TRANSACTIONS  A 


43 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Fig.  4 — Evolution  of  the  (a)  bulk  PDAS  and  ( b )  standard  deviation  of  the  bulk  PDAS  as  a  function  of  the  length  scale.  The  error  bars  are  for 
20  different  microstructure  instantiations  for  cubic,  hexagonal,  and  random  patterns  with  a  noise  fraction  of  0.40ao-  The  VWR  technique  (a  =  2) 
and  the  Voronoi  technique  (dCT jt  =  0.12)  were  used  to  calculate  the  bulk  PDAS  values  for  each  microstructure. 


pattern  type  and  the  measurement  technique.  With  both 
measurement  techniques,  the  cubic  pattern  is  closest  to 
the  expected  PDAS  (without  noise),  followed  by  the 
hexagonal  pattern  and  last  the  random  pattern.  For  all 
pattern  types,  the  Voronoi  technique  was  closer  to  the 
expected  PDAS  than  the  VWR  method.  Based  on  these 
tests,  the  80  x  80  patterns  are  chosen  to  illustrate  the 
subsequent  trends  since  multiple  microstructure  instan¬ 
tiations  should  not  be  needed  to  examine  the  subsequent 
trends  in  terms  of  microstructure  and  local  PDAS 
measurement  technique. 

B.  Bulk  PDAS  and  Coordination  Number  Sensitivity 
Analysis 

The  bulk  PDAS  and  bulk  coordination  number  was 
calculated  for  multiple  microstructures  to  understand 
the  parameter  sensitivity  for  the  two  different  local 
PDAS  measurement  techniques.  The  a  parameter  was 
varied  for  the  VWR  method  (Figure  5)  and  the  dCI ;t 
parameter  was  varied  for  the  Voronoi  method  (Fig¬ 
ure  6).  The  expected  values  are  represented  by  the 
dotted  lines  in  all  subfigures  (except  Figure  5(d)  where 
the  bulk  bulk  coordination  numbers  significantly  un¬ 
derestimate  the  expected  coordination  number). 

There  are  several  trends  that  are  evident  from 
Figures  5  and  6.  In  the  case  of  the  VWR  method,  both 
the  bulk  PDAS  and  the  bulk  coordination  number  are 
affected  by  the  parameter  at,  as  has  been  noted  before. 
As  at  is  increased,  the  bulk  PDAS  value  increases 
towards  the  expected  PDAS  of  1  and  the  bulk  coordi¬ 
nation  number  increases  for  all  microstructures  exam¬ 
ined  in  this  study.  As  the  noise  level  is  increased  for  both 
the  cubic  and  hexagonal  microstructures,  the  trend  is  for 
the  bulk  PDAS  to  deviate  further  from  the  expected 
value.  However,  for  high  values  of  tz,  the  bulk  PDAS  for 
cubic  microstructures  with  higher  noise  levels  does  begin 

METALLURGICAL  AND  MATERIALS  TRANSACTIONS  A 


to  approach  the  expected  PDAS  at  a  faster  rate  than 
that  of  the  microstructures  with  lower  noise  levels 
(Figure  5(a)).  This  behavior  is  a  function  of  two 
interacting  effects:  (1)  higher  a  values  admitting  more 
nearest  neighbors,  which  have  larger  spacings  (as  can  be 
observed  in  Figure  5(c)),  and  (2)  higher  noise  levels 
increase  the  standard  deviation  of  the  nearest  neighbor 
spacings,  which  (in  combination  with  a)  can  result  in 
admitting  even  more  nearest  neighbors.  Interestingly, 
the  bulk  coordination  number  for  the  hexagonal  mi¬ 
crostructures  in  Figure  5(d)  does  not  approach  the 
expected  value  of  6  for  any  of  the  a  values  chosen.  This 
behavior  helps  to  explain  the  larger  deviation  of  this 
technique  from  the  expected  PDAS  value  for  hexagonal 
microstructures  in  Figure  5(b).  However,  the  inability  of 
this  technique  to  capture  the  bulk  coordination  number 
even  for  small  levels  of  random  noise  in  hexagonal 
microstructures  (0.05ao)  signifies  that  while  it  may  be 
useful  for  exploring  some  trends  between  different 
microstructures,  it  may  significantly  under-predict  co¬ 
ordination  numbers  for  hexagonal  dendritic  configura¬ 
tions  (which  may  skew  the  PDAS  statistics  as  well). 

In  the  case  of  the  Voronoi  method,  both  the  bulk 
PDAS  and  the  bulk  coordination  number  are  similarly 
affected  by  the  parameter  r/crjt.  As  dCI is  increased,  the 
trend  is  for  both  the  bulk  PDAS  and  bulk  coordination 
number  to  either  decrease  or  remain  constant.  The  noise 
level  of  the  cubic  or  hexagonal  microstructure  sig¬ 
nificantly  affects  the  behavior  of  these  responses  with 
respect  to  da jt.  For  low  noise  levels,  the  Voronoi  method 
with  a  large  da jt  value  (dCI jt  >  0.04)  captures  both  the 
expected  PDAS  value  and  the  expected  coordination 
number — this  behavior  indicates  that  this  method  is 
capable  of  capturing  the  expected  behavior  as  the  noise 
is  reduced  to  zero,  in  contrast  to  the  VWR  method.  The 
only  unusual  behavior  at  low  noise  levels  for  the 
Voronoi  method  is  that  obtained  for  small  da jt  values 
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Fig.  5 — Evolution  of  the  (a,  b )  bulk  PDAS  and  (c,  d)  bulk  coordination  number  as  a  function  of  the  parameter  a  for  the  Voronoi  Warnkens- 
Reed  method  for  both  (a,  c)  cubic  and  (b,  d)  hexagonal  microstructures  (80  x  80)  with  different  noise  levels.  The  dotted  line  represents  the  ex¬ 
pected  value  for  a  zero-noise  microstructure. 


(rfCrit  <  0.02).  For  rfcri,  =  0,  the  bulk  coordination  num¬ 
ber  for  the  cubic  microstructure  is  6,  regardless  of  the 
noise  level,  because  the  4  nearest  neighbors  always  form 
a  face  on  the  Voronoi  polygon,  while,  on  average,  2  of 
the  4  corner  neighbors  form  another  face  of  the  polygon 
(since  there  is  no  exclusion  of  any  of  the  faces  (i.e.,  the 
dCI ^  parameter),  the  average  coordination  number  is 
always  6).  In  summary,  the  ability  to  capture  the 
expected  behavior  (i.e.,  zero-noise,  ideal  microstruc¬ 
tures)  in  microstructures  with  small  noise  levels  provides 
some  degree  of  confidence  that  the  Voronoi  method 
with  r/crit  >  0.04  is  a  technique  that  may  be  appropriate 
for  examining  local  PDAS  and  coordination  number. 
Unfortunately,  the  bulk  coordination  number  (and, 
consequently,  the  bulk  PDAS)  still  changes  as  a  function 
of  noise  in  the  microstructure  (Figures  6(c)  and  (d))  and 
there  is  no  “true”  coordination  number  to  evaluate  what 


exact  value  of  r/crjt  is  appropriate.  In  practice,  it  should  be 
realized  that  the  choice  of  the  dCT ;t  value  with  the  Voronoi 
technique  can  affect  both  the  bulk  PDAS  behavior  and 
the  bulk  coordination  number,  but  should  be  appropriate 
for  analyzing  trends  in  microstructures.  In  the  subsequent 
subsections,  the  local  PDAS  and  CN  calculations  over  a 
number  of  different  <fcrjt  values  for  the  Voronoi  technique 
will  be  explored  to  understand  their  correlations  and  the 
evolution  of  these  distributions. 

C.  Local  PDAS  and  Coordination  Number  Sensitivity 
Analysis 

The  local  PDAS  and  coordination  number  distribu¬ 
tions  cause  the  observed  changes  in  the  bulk  behavior. 
To  better  understand  how  these  distributions  change  as 
a  function  of  the  r/crj,  value  for  the  Voronoi  technique 
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Fig.  6 — Evolution  of  the  (a,  b )  bulk  PDAS  and  (c,  d)  bulk  coordination  number  as  a  function  of  the  parameter  dcr jt  for  the  Voronoi  method  for 
both  (a,  c)  cubic  and  (b,  d)  hexagonal  microstructures  (80  x  80)  with  different  noise  levels.  The  dotted  line  represents  the  expected  value  for  a 
zero-noise  microstructure. 


(and  for  a  with  the  VWR  technique,  not  shown),  six 
different  noise  levels  for  cubic  and  hexagonal  80  x  80 
microstructures  were  generated  and  seven  different  dcr jt 
values  were  performed  for  each  microstructure.  As  an 
example  of  the  probability  distribution  plots  used  to 
show  the  evolution  of  local  PDAS  and  coordination 
with  noise.  Figure  7  shows  their  evolution  with  dcr jt 
(multiple  lines  within  each  subfigure)  for  a  random 
microstructure.  There  are  several  features  within  each 
plot  meant  to  aid  in  their  analysis  and  comparison. 
First,  each  plot  displays  the  local  PDAS  probability 
distribution  (left  plot,  red  lines)  and  the  local  coordina¬ 
tion  number  histogram  (right  plot,  red  bars).  To 
delineate  the  different  values  of  rfcrit ,  i/crit  =  0  is  shown 
as  a  thick  maroon  line  (or  thick  maroon  bar  for  the 
coordination  number  plot)  and  the  remainder  of  the 


lines  (bars)  indicate  monotonic  increases  in  increments 
of  rfcrit  by  0.02.  On  each  plot,  the  bulk  PDAS/coordi- 
nation  numbers  for  the  Voronoi  method  with  various 
r/crit  values  are  shown  as  blue  lines,  which  are  also 
incremented  from  0.00  (right,  thick  line)  to  0.12  (left)  in 
increments  of  0.02.  Since  in  some  cases  these  lines  may 
overlap,  a  circle  marker  symbol  is  attached  to  the  line, 
where  the  size  of  the  symbol  corresponds  to  the 
magnitude  of  £/crjt.  For  the  case  of  the  random  mi¬ 
crostructure,  increases  in  dCI ;t  decrease  the  local  and  bulk 
PDAS  distribution  as  well  as  the  coordination  number 
distribution.  Furthermore,  as  ordered  microstructures 
( e.g .,  cubic  and  hexagonal)  approach  a  more  disordered 
state,  the  coordination  number  distribution  is  more 
likely  to  broaden  and  assume  a  distribution  of  different 
local  coordination  numbers  (with  4,  5,  and  6  being  the 
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Local  PDAS  Local  Coordination  Number 


Fig.  7 — Cumulative  probability  plot  of  the  local  PDAS  (left  plot)  and  coordination  number  (right  plot)  distribution  as  a  function  of  the  pa¬ 
rameter  r/ent  for  the  Voronoi  method  for  a  random  microstructures  (80  x  80).  The  red  (dotted  blue)  line  represents  the  local  PDAS/coordination 
number  distribution  (bulk  PDAS/coordination  number)  for  dCI a  values  of  0.00  (right,  thick  line)  to  0.12  (left)  in  increments  of  0.02. 


most  frequent).  The  remainder  of  this  section  will  focus 
on  the  evolution  of  the  local  statistics  for  cubic  and 
hexagonal  microstructures  with  varying  degrees  of 
noise. 

1 .  Voronoi  technique  with  cubic  microstructures 

The  local  PDAS  and  coordination  number  distribu¬ 
tions  cause  the  observed  changes  in  the  bulk  behavior. 
To  better  understand  how  these  distributions  change  as 
a  function  of  the  da\t  value  for  the  Voronoi  technique 
(and  for  a  with  the  VWR  technique,  not  shown),  six 
different  noise  levels  for  cubic  and  hexagonal  80  x  80 
microstructures  were  generated  and  seven  different  da ;t 
values  were  performed  for  each  microstructure.  Figure  8 
shows  the  evolution  of  the  local  PDAS  and  coordination 
number  distributions  as  a  function  of  noise  level 
(Figures  8(a)  through  (f))  and  da ;t  (multiple  lines  within 
each  subfigure).  Each  subfigure  displays  the  local 
PDAS/coordination  number  statistics  in  the  same  man¬ 
ner  as  that  described  for  the  random  microstructure 
(Figure  7).  The  hypothesis  for  the  following  set  of  cubic 
microstructures  is  that  as  the  noise  is  increased,  the 
standard  deviation  of  the  local  PDAS/coordination 
number  distributions  ( i.e .,  the  slope  of  the  distribution 
functions)  increases,  but  the  mean  values  will  tend  to 
reflect  that  of  the  base  cubic  microstructure  (i.e., 
PDAS  =  1  and  bulk  coordination  number  of  4) — the 
values  of  dCI\ t  that  best  capture  this  might  best  describe 
the  local  statistics. 

Several  trends  are  noticeable  from  the  various 
probability  distribution  function  plots  in  Figure  8.  The 
first  observation  is  that  at  low  noise  levels  (e.g.,  0.05«o) 
and  low  dCI jt  values,  the  PDAS  deviates  greatly  from  the 
standard  S-shaped  normal  distribution  curve.  This 
abnormal  distribution  can  be  further  explained  by  the 
local  coordination  number  distribution  (right  plots, 
Figure  8).  The  local  coordination  number  plots  in 
Figure  8  indicate  that  these  conditions  cause  a  sig¬ 
nificant  number  of  local  coordination  numbers  above  4 
to  appear  (bulk  coordination  number  as  high  as  6,  as 
noted  previously),  even  though  the  low  noise  levels 
should  approach  coordination  numbers  similar  to  a 
zero-noise  cubic  microstructure  (i.e.,  bulk  coordination 
number  of  4).  Since  there  is  very  little  noise,  these 
nearest  neighbors  are  likely  the  diagonal  second-nearest 
neighbors,  which  have  a  higher  spacing  than  the  four 


nearest  neighbors  in  the  cubic  arrangement;  thus,  the 
local  PDAS  values  are  heavily  influenced  by  the 
inclusion  of  the  5th,  6th,  and  perhaps  higher-numbered 
neighbors. 

A  second  observation  is  that  the  distributions  become 
more  similar  in  shape  and  symmetric  about  the  mean  as 
the  noise  is  increased,  which  reflects  that  the  standard 
deviations  of  the  distributions  are  similar  and  the 
distributions  are  merely  offset  by  different  mean  values, 
which  linearly  correlate  to  dCI\t  at  higher  noise  levels.  It 
should  be  noted  that,  for  every  noise  level,  the  higher 
values  of  dCI ;t  (={0.10,0.12})  tend  to  approach 
PDAS  =  1 — the  PDAS  of  the  zero-noise  cubic  mi¬ 
crostructure.  The  slope  of  the  local  PDAS  probability 
distribution  is  also  decreased  with  higher  noise  levels  for 
dent  (={0.10,0.12}),  i.e.,  increased  noise,  increased 
standard  deviation  of  the  local  PDAS  distribution. 

A  third  observation  is  that,  as  the  noise  level  is 
increased,  more  local  coordination  numbers  lower  than 
4  are  admitted.  This  statement  is  particularly  true  for 
higher  values  of  dCI jt;  hence,  as  t/crjt  is  increased,  there  is  a 
higher  probability  of  admitting  local  coordination 
numbers  of  either  2  or  3.  Conversely,  as  dCI ;t  is 
decreased,  there  is  a  higher  probability  of  admitting 
local  coordination  numbers  of  6  or  greater.  In  terms  of 
the  mean  behavior,  this  results  in  bulk  coordination 
numbers  of  4  and  6  for  z/cHt  =  0,12  and  z/crit  =  0.00, 
respectively.  So,  for  cubic  microstructures  with  varying 
degrees  of  noise,  z/crit  values  of  0.10  or  0.12  yields  local 
PDAS  and  coordination  numbers  that  are  consistent 
with  the  zero-noise  cubic  microstructure,  while  display¬ 
ing  an  increase  in  the  standard  deviation  of  the 
distribution  with  increasing  noise,  as  would  be  expected. 

2.  Voronoi  technique  with  hexagonal  microstructures 

The  same  sort  of  probability  distribution  plots  are 
shown  for  hexagonal  microstructures  in  Figure  9. 
Again,  the  hypothesis  is  that  as  the  noise  within  the 
microstructure  is  increased,  the  PDAS  should  still  trend 
towards  1  and  the  local  coordination  number  should  be 
approximately  6,  but  the  standard  deviation  of  the 
distribution  will  increase.  There  are  a  few  interesting 
observations  that  contrast  those  with  the  cubic  mi¬ 
crostructures. 

First,  in  the  case  of  low  values  of  <fcrjt,  there  is  not  the 
same  abnormal  behavior  at  low  noise  levels  as  with  the 
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Fig.  8 — Cumulative  probability  plot  of  the  local  PDAS  (left  plot)  and  coordination  number  (right  plot)  distribution  as  a  function  of  the  pa¬ 
rameter  dCT it  for  the  Voronoi  method  for  cubic  microstructures  (80  x  80)  with  different  noise  levels:  ( a )  0.05ao,  ( b )  O.lOflo,  (c)  0.20ao,  (d)  0.30ao, 
( e )  0.40tfo,  and  if)  0.50tfo-  The  red  (dotted  blue)  line  represents  the  local  PDAS/coordination  number  distribution  (bulk  PDAS/coordination  num¬ 
ber)  for  *4^  values  of  0.00  (right,  thick  line)  to  0.12  (left)  in  increments  of  0.02. 
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Fig.  9 — Cumulative  probability  plot  of  the  local  PDAS  (left  plot)  and  coordination  number  (right  plot)  distribution  as  a  function  of  the  pa¬ 
rameter  dfcrjt  for  the  Voronoi  method  for  hexagonal  microstructures  (80  x  80)  with  different  noise  levels:  ( a )  0.05«o,  ( b )  O.lOflo,  (c)  0.20ao,  (d) 
0.30ao,  (e)  0.40tfo,  and  (/)  0.50«o-  The  red  (dotted  blue)  line  represents  the  local  PDAS/coordination  number  distribution  (bulk  PDAS/coordina- 
tion  number)  for  dcr jt  values  of  0.00  (right,  thick  line)  to  0.12  (left)  in  increments  of  0.02. 
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cubic  microstructures.  This  is  perhaps  an  expected 
result,  because  the  six  nearest  neighbors  effectively 
shield  the  second-nearest  neighbors,  which  enabled  the 
abnormal  behavior  in  the  distributions  for  the  cubic 
microstructures.  For  the  same  small  noise  levels  with  the 
hexagonal  microstructures,  there  is  also  very  little 
difference  between  the  local  PDAS  or  the  local  coordi¬ 
nation  number  distributions  compared  to  that  observed 
within  the  cubic  microstructures. 

Second,  as  the  noise  in  increased  within  the  hexagonal 
microstructure,  the  spread  of  the  bulk  PDAS  and  bulk 
coordination  numbers  increased.  However,  in  contrast 
to  the  cubic  microstructures,  the  bulk  PDAS  values  drop 
below  PDAS  =  1  and  a  lower  value  of  rfcrjt  =  0.04  or 
^crit  =  0.06  is  actually  closer  to  the  ideal-hexagonal- 
microstructure  PDAS  value. 

Third,  the  tendency  is  for  lower  coordination  num¬ 
bers  to  be  calculated  with  increasing  noise  in  the 
microstructure  and  with  increasing  «fCrit  values.  For  the 
highest  noise  level  explored  herein,  the  extreme  bulk 
coordination  numbers  are  4  and  6  for  da ;t  =  0.12  and 
^crit  =  0.00,  respectively.  Interestingly,  these  bulk  coor¬ 
dination  numbers  are  nearly  identical  to  those  of  the 
cubic  microstructures  with  the  highest  noise  level — so 
there  may  be  a  convergence  in  the  behavior  of  the 
Voronoi  technique  with  respect  to  the  higher  noise 
levels  of  the  microstructure.  For  smaller  noise  levels  in 
the  cubic  and  hexagonal  microstructures,  dCI ;t  values 
between  0.04  and  0.12  give  reasonable  behavior  in 
comparison  to  the  known  values  from  a  zero-noise 
microstructure. 

3.  Voronoi  Warnken—Reed  technique 

The  same  analysis  was  also  conducted  to  examine 
the  variability  due  to  the  noise  level  and  the  parameter 
a  within  the  VWR  technique.  The  subsequent  plots  are 
attached  in  Appendix.  What  is  immediately  apparent  is 
that  a  large  fraction  of  points  have  a  calculated 
coordination  number  of  3  using  this  technique,  irre¬ 
spective  of  noise  level  and  even  a  value.  Recall  that  this 
is  an  initial  parameter  that  is  set  within  the  VWR 
method,  whereby  the  three  nearest  neighbors  were 
initially  selected  to  calculate  the  standard  deviation  and 
then  this  value  is  used  with  a  to  decide  on  the  next- 
closest  neighbor — an  initial  value  of  2  is  too  few  for  a 
standard  deviation  calculation  and  a  value  of  4 
eliminates  the  possibility  of  local  coordination  numbers 
with  values  less  than  4.  Moreover,  the  bulk  coordina¬ 
tion  number  never  exceeds  a  value  of  4.5,  even  with 
a  =  2.0,  for  hexagonal  structures  (even  at  low  noise 
levels).  This  finding  indicates  that  this  distance-based 
technique  is  certainly  underestimating  the  coordination 
number  in  hexagonal  microstructures,  i.e.,  compare 
Figure  A2  in  Appendix  with  Figure  9  to  compare 
techniques  for  the  exact  same  microstructure.  In  a 
distance-based  technique  such  as  this,  the  addition  of 
these  nearest  neighbors  always  increases  the  local 
PDAS.  Hence,  the  larger  a  is  for  the  VWR,  the  more 
skewed  the  local  PDAS  distribution  is  towards  higher 
local  PDAS  values. 
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D.  Local  Correlation  for  Different  Parameters 

The  subsequent  datasets  can  also  be  analyzed  to 
examine  the  correlation  between  local  PDAS  and  local 
coordination  numbers.  This  sort  of  analysis  can  help 
answer  the  question  “is  there  a  high  degree  of  correla¬ 
tion  between  the  local  PDAS  values  for  different 
technique  parameter  values?”  Another  valid  question 
is  that  of  "is  there  a  high  degree  of  correlation  between 
the  local  PDAS  values  and  the  local  coordination 
numbers  for  different  technique  parameter  values?” 
Both  of  the  aforementioned  questions  are  subsequently 
addressed  in  the  next  two  subsubsections. 

1 .  Local  PDAS  correlation 

Figure  10  shows  a  few  examples  of  correlations 
between  the  local  PDAS  values  for  two  different 
methods — (a)  Voronoi  method  and  (b)  VWR  method — 
using  two  different  parameter  values  (r/cr;t  of  {0.00, 0. 12} 
and  a  of  {1.0, 2.0},  respectively).  These  scatter  plots 
combine  all  of  the  various  80  x  80  cubic  and  hexagonal 
microstructures  explored  in  Section  III-C.  The  scatter 
plots  were  incrementally  plotted  such  that  the  largest 
noise  level  was  plotted,  then  the  second  largest,  etc.  In 
this  manner,  the  lower-noise  data  points  are  plotted  over 
higher-noise  data  points,  but  overall  this  plot  identifies 
the  degree  of  randomness  or  bias  between  the  two 
parameters  chosen  for  the  same  technique.  Note  that 
each  data  point  represents  the  technique  being  applied 
to  the  same  point  with  the  same  local  environment. 

There  are  some  similarities  and  differences  for  the 
correlation  plots  in  Figure  10.  Since  the  colors  in  these 
plots  indicate  the  different  noise  levels,  it  is  observed 
that  the  larger  noise  levels  result  in  a  larger  spread  of  the 
PDAS  values,  as  would  be  expected  from  the  previous 
analyses  conducted  herein.  In  both  cases,  though,  there 
is  not  a  linear  correlation  between  the  two  different 
parameter  levels,  i.e..  knowing  the  local  PDAS  at  one 
point  for  a  given  parameter  level  is  insufficient  infor¬ 
mation  to  predict  the  local  PDAS  at  that  same  point 
with  a  different  parameter  level  for  the  same  technique. 

Some  of  the  differences  are  as  follows.  First,  notice 
that  the  Voronoi  technique  in  Figure  10(a)  has  some 
values  that  lie  above  the  45  deg  line.  This  indicates  that 
in  some  cases,  the  Voronoi  technique  with  a  higher  r/crj, 
will,  in  fact,  exclude  some  closer  local  PDAS  values, 
which  causes  the  local  PDAS  value  to  increase  with 
increasing  rfCnt-  However,  this  situation  is  not  a  common 
occurrence,  as  can  be  observed  from  the  majority  of 
points  that  lie  below  this  line.  This  behavior  contrasts 
the  VWR  technique  in  Figure  10,  which  has  all  local 
PDAS  values  either  on  the  45  deg  line  or  above  it.  Of 
course,  this  is  as  expected  since  the  VWR  is  a  distance- 
based  technique,  so  an  increased  a  results  in  either  the 
same  number  of  nearest  neighbors  (and  local  PDAS)  or 
a  higher  number  of  neighbors  (with  a  larger  local 
PDAS).  Second,  at  low  noises  (blue  and  light  blue  data 
points),  the  Voronoi  method  starts  to  deviate  from  the 
expected  linear  correlation  between  the  two  clan  values. 
This  behavior  is  caused  by  dCI jt  =  0.0  counting  nearest 
neighbors  that  have  moved  further  away,  thus  biasing 
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Fig.  10 — Scatter  plot  of  the  correlation  between  the  local  PDAS  for  two  different  methods — (a)  Voronoi  method  and  ( b )  Voronoi  Warnken- 
Reed  method — using  two  different  parameter  values  (dCT jt  of  {0.00,0.12}  and  a  of  (1.0, 2.0},  respectively).  These  scatter  plots  combine  all  of  the 
various  80  x  80  cubic  and  hexagonal  microstructures  explored  in  Section  III— C  and  the  colors  correspond  to  the  different  microstructure  noise 
levels  (see  legend  on  right). 
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Fig.  11 — Scatter  plot  of  the  correlation  between  the  local  PDAS  values  and  the  local  coordination  numbers  for  Voronoi  method  with  two  differ¬ 
ent  parameter  values:  dCTn  of  (a)  0.04  and  ( b )  0.12.  These  scatter  plots  combine  all  of  the  various  80  x  80  cubic  and  hexagonal  microstructures 
explored  in  Section  III— C  and  the  colors  correspond  to  the  different  microstructure  noise  levels  (see  legend  on  right).  The  larger  offset  dots  indi¬ 
cate  the  bulk  PDAS  values  as  a  function  of  noise  level  and  the  (nearest)  coordination  number. 


the  local  PDAS  by  including  these  points — this  is  espe¬ 
cially  evident  for  the  case  of  cubic  microstructures  where 
second-nearest  neighbors  are  counted  for  dCT it  =  0.0.  On 
the  other  hand,  the  higher  r/crjt  value  instead  excludes  these 
further  points  and  maintains  a  local  PDAS  of  1.0,  in 
agreement  with  the  zero-noise  microstructure. 

2.  Local  PDAS  and  local  coordination  number 
correlation 

Figure  1 1  shows  a  few  examples  of  correlations 
between  the  local  PDAS  values  and  the  local  coordina¬ 
tion  numbers  for  Voronoi  method  with  two  different 
parameter  values  (dCI jt  of  {0.04,0.12}).  These  scatter 
plots  are  similar  to  Figure  10  in  that  they  combine  all  of 
the  various  80  x  80  cubic  and  hexagonal  microstruc¬ 


tures  explored  in  Section  III  C.  The  scatter  plots  were 
incrementally  plotted  form  the  largest  noise  level  to  the 
smallest  in  a  similar  manner.  However,  the  coordination 
numbers  are  discrete  values,  so  the  scatter  plot  takes  on 
a  different  form.  In  order  to  understand  how  the  bulk 
PDAS  is  also  evolving  within  each  coordination  num¬ 
ber,  a  larger  circle  marker  symbol  is  included  above  each 
discrete  coordination  number  to  indicate  this  trend. 

There  are  some  trends  that  emerge  from  the  correla¬ 
tion  plots  in  Figure  11.  In  general,  increasing  local 
coordination  number  results  in  higher  local  and  bulk 
PDAS  values.  Local  coordination  numbers  of  4,  5,  and  6 
are  common  over  all  noise  levels  and  both  dCI j,  values, 
with  some  lower  and  higher  local  coordination  numbers 
occurring  with  an  increase  in  the  noise  of  the  mi- 
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crostructure.  In  general,  the  higher  dalt  values  results  in 
lower  local  coordination  numbers  and,  consequently, 
lower  PDAS  values.  In  terms  of  the  bulk  PDAS 
behavior  with  respect  to  the  individual  coordination 
numbers,  two  coordination  number  regimes  are  ob¬ 
served:  (1)  For  coordination  numbers  less  than  5,  the 
bulk  PDAS  decreases  with  increasing  noise  level,  and 
(2)  for  coordination  numbers  greater  than  5,  the  bulk 
PDAS  increases  with  increasing  noise  level.  In  both 
cases,  the  general  trend  is  that  the  bulk  PDAS  for  each 
coordination  number  will  move  further  away  from 
PDAS  =  1 .0  with  increasing  noise  in  the  microstruc¬ 
ture.  The  only  exception  is  for  the  coordination 
number  of  5,  which  is  a  function  of  the  competing 
influences  from  the  cubic  (coordination  number 
of  4)  and  hexagonal  (coordination  number  of  6) 
microstructures. 


IV.  CONCLUSIONS 

In  this  work,  four  series  of  simulations  using  cubic 
and  hexagonal  synthetic  microstructures  with  varying 
degrees  of  noise  were  utilized  to  study  the  effects  of 
length  scale,  base  microstructure,  microstructure  vari¬ 
ability,  and  technique  parameters  on  the  local  PDAS 
distribution,  local  coordination  number  distribution, 
bulk  PDAS,  and  bulk  coordination  number  (calculated 
as  the  mean  of  the  local  statistics).  The  length  scale 
study  used  20  instantiations  to  identify  a  size  of  our 
synthetic  microstructure  whereby  the  variability  in  the 
mean  statistics  starts  to  converge — an  80  x  80  mi¬ 
crostructure  (i.e.,  6400  dendrites)  was  chosen  herein. 
While  a  number  of  the  quantitative  findings  are  detailed 
herein,  a  few  important  findings  help  to  evaluate  the 
relative  capability  of  these  approaches  to  quantify  theses 
local  statistics. 

1.  For  the  Voronoi  tesselation  technique,  utilizing  the 
tfcrit  (polygon-side-length)  criterion  (dcrjt  >  0.04)  was 
necessary  to  correctly  capture  bulk  coordination 
numbers  for  cubic  and  hexagonal  microstructures 
with  small  amounts  of  disorder  (noise). 

2.  For  the  VWR  method,  the  a  parameter  in  the  range 
examined  herein  is  not  able  to  capture  the  bulk 
coordination  numbers  for  hexagonal  microstruc¬ 
tures  with  small  amounts  of  noise,  signifying  that  it 
may  significantly  under-predict  coordination  num¬ 
bers  for  hexagonal  dendritic  configurations  (which 
results  in  under-predicting  the  PDAS  statistics  as 
well). 

3.  The  subsequent  analysis  of  the  local  coordination 
numbers  and  local  PDAS  distributions  identifies  a 
number  of  trends  in  the  local  distributions  that  help 
explain  the  behavior  of  these  techniques  and  their 
bulk-calculated  statistics.  In  general,  the  right 


trends  are  observed — with  increasing  noise  in  the 
microstructures,  the  Voronoi  technique  correctly 
captures  the  increasing  standard  deviation  of  the  lo¬ 
cal  PDAS  distribution. 

4.  In  terms  of  the  correlation  study,  there  is  not  a 
clear  linear  correlation  between  different  da jt  values 
for  the  local  statistics;  this  may  be  partly  associated 
with  the  general  (weak,  but  clearly  evident)  trend 
that  an  increase  in  the  local  coordination  number 
tends  to  increase  the  local  PDAS  value. 

Overall,  the  Voronoi  tesselation  technique  with  the 
dcrit  (polygon-side-length)  criterion  best  captures  the 
wide  range  of  directionally  solidified  dendritic  mi¬ 
crostructures.  This  systematic  study  of  the  different 
techniques  for  quantifying  local  primary  dendrite  arm 
spacings  evaluates  their  capability  to  capture  this 
important  microstructure  feature  in  different  syn¬ 
thetically  derived  dendritic  microstructures,  which 
can  be  an  important  step  for  correlating  with  both 
processing  and  properties  in  directionally  solidified 
microstructures,  such  as  in  single  crystal  nickel-based 
superalloys. 
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APPENDIX 

Figures  A1  and  A2  present  data  pertaining  to  how  the 
local  PDAS  and  coordination  number  distributions 
change  as  a  function  of  the  a  value  for  the  VWR 
technique.  Six  different  noise  levels  for  cubic  and 
hexagonal  80  x  80  microstructures  were  generated 
(Figures  A 1  and  A2,  respectively)  and  11  different  a 
values  were  performed  for  each  microstructure.  The  a 
values  were  incremented  from  1.0  to  2.0  in  increments  of 
0.1  for  the  subsequent  analysis.  All  other  aspects  are 
similar  in  nature  to  Figures  7,  8,  and  9,  which  show  the 
evolution  of  these  local  distributions  with  dCI for  the 
Voronoi  technique.  Further  discussion  of  these  plots 
and  their  results  is  in  Section  III  C  3. 
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Fig.  A1 — Cumulative  probability  plot  of  the  local  PDAS  (left  plot)  and  coordination  number  (right  plot)  distribution  as  a  function  of  the  pa¬ 
rameter  a  for  the  Voronoi  Warnken-Reed  method  for  cubic  microstructures  (80  x  80)  with  different  noise  levels:  ( a )  0.05ao,  ( b )  O.lOflo,  (c) 
0.20ao,  (d)  0.30«o,  (e)  0.40ao,  and  (/)  0.50flo-  The  red  (dotted  blue)  line  represents  the  local  PDAS/coordination  number  distribution  (bulk 
PDAS/coordination  number)  for  a  values  of  0.00  (right,  thick  line)  to  0.12  (left)  in  increments  of  0.02. 
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Fig.  A2 — Cumulative  probability  plot  of  the  local  PDAS  (left  plot)  and  coordination  number  (right  plot)  distribution  as  a  function  of  the  pa¬ 
rameter  a  for  the  Voronoi  Warnken-Reed  method  for  hexagonal  microstructures  (80  x  80)  with  different  noise  levels:  (a)  0.05ao,  ( b )  O.lOtro, 
(c)  0.20flo,  (d)  0.30#o,  (e)  0.40flo,  and  (/)  0.50«o-  The  red  (dotted  blue)  line  represents  the  local  PDAS/coordination  number  distribution  (bulk 
PD AS/coordination  number)  for  a  values  of  0.00  (right,  thick  line)  to  0.12  (left)  in  increments  of  0.02. 
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